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Abstract 

By a combination of analytical and numerical techniques, we analyze the continuum limit of the 
integrable 3(g)3®3®3. . . sl{2/l) superspin chain. We discover profoundly new features, including 
a continuous spectrum of conformal weights, whose numerical evidence is infinite degeneracies of 
the scaled gaps in the thermodynamic limit. This indicates that the corresponding conformal 
field theory has a non compact target space (even though our lattice model involves only finite 
dimensional representations). We argue that our results are compatible with this theory being the 
level k = 1, 'S(7(2/l) WZW model' (whose precise definition requires some care). In doing so, 
we establish several new results for this model. With regard to potential applications to the spin 
quantum Hall efi'ect, we conclude that the continuum limit of the 3Cg)3® 3®3 . . . sZ(2/l) integrable 
superspin chain is not the same as (and is in fact very different from) the continuum limit of the 
corresponding chain with two-superspin interactions only, which is known to be a model for the 
spin quantum Hall effect. The study of possible RG flows between the two theories is left for 
further study. 

1 Introduction 

The sigma model approach [1] to phase transitions in non interacting disordered systems provides a 
convenient and appeahng description of the physics at hand. The best known example is the tran- 
sition between plateaux in the integer quantum Hall effect (class A), which can be described using 
a U{2n)/U{n) x U{n) sigma model in the limit n ^ 0. More recently, the equivalent 'spin quan- 
tum Hall effect' in d-wave superconductors has been considered (class C [2]), and described with the 
SP{2n)/U{n) sigma model, in the limit n 0[3]. 

In both cases, the existence of a delocalization transition is associated with the masslessness in the 
IR of the sigma model at topological angle 6 = n. Identification of the corresponding conformal field 
theory would lead to a presumably exact determination of critical exponents which have been studied 
in great details in numerical works [4], [5] as well as in experiments. 

Unfortunately, the determination of the IR fixed points in these problems has proven extremely 
difficult. The canonical example of what 'might' happen is the case n = 1 of either replica model, 
which coincides with the 0(3) sigma model at 6* = tt. In this case, there is a wealth of evidence [6] that 
the IR fixed point is the SU{2) Wess Zumino model at level A; = 1. In the fiow, the field is promoted 
from being coset valued to group valued, a highly non perturbative feature which might have been 
hard to believe, were it not for the existence of exact solutions [7]. 

The most reasonable way to understand what happens in the disordered system is first to trade 
the replica hmit for a super group [8], [9], U{1, l/2)/J7(l/l) x U{1/1) in the quantum Hall effect, and 
OSP{2n/2n)/U{n\n) in the spin quantum Hall effect. Once the corresponding sigma models have been 
satisfactorily defined [10], it is then tempting to conjecture, by analogy with the 0(3) example above, 
what the IR fixed point might be. In this way, Zirnbauer [11] arrived at a proposal for the quantum Hall 
effect that is related with the WZW model (at level fc = 1) on PSU{1, 1/2). Roughly, this is obtained 
by observing that the base of the supersymmetric target space is x S"^, and then guessing that this 
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gets 'promoted' to x in the IR. One then tries to build the minimal proposal in agreement with 
this and known results on the exponents. It is important that for the PSU{1, 1/2) WZW model the 
kinetic term is an exactly marginal deformation [12], [13]. The final proposal of Zirnbaucr corresponds 
to a particular value of this kinetic term, and is not the WZW model. Hence, the symmetry is a global 
symmetry Gl x Gr instead of a local one. Unfortunately it is not clear whether this proposal gives 
the correct IR fixed point. In particular, it is rather unpleasant to have to select a particular value 
for an exactly marginal coupling constant if one hopes to describe what seems to be genuine universal 
physics. 

A less abstract direction of attack comes from network models. These can be shown to correspond, 
in the proper anisotropic limit to quantum spin chains with super group symmetries [14]. Ideally, 
the problem would then be solved if one were able to apply in this situation the formidable arsenal 
developed over the years in the study of quantum spin chains. Unfortunately, there are difficulties here 
as well. For the ordinary as well as the spin quantum Hall effect, the Hilbert space of the chain is of 
type VCSiV^V^V^---, while the Hamiltonian is the invariant Casimir, the proper generalization of 
S.S. In the ordinary quantum Hall effect, V and V arc infinite dimensional, while in the spin quantum 
Hall effect, V and V are the two conjugate three dimensional representations. 

In either case, the Hamiltonian is not integrable. More precisely, the Hamiltonian does not seem to 
derive from a family of commuting transfer matrices obtained via a standard quantum group approach. 
In the spin quantum Hall case, the Hamiltonian can however be diagonalized [15] (more precisely, the 
zeroes of the characteristic polynomial can be found, as the Hamiltonian itself is not fully diagonal- 
izable) using known results about the XXZ chain and representation theory of the Temperley Lieb 
algebra [16]. No such approach seems possible in the ordinary quantum Hall spin chain. 

Since one is after the universality class of the IR fixed point, it is natural to wonder [17] whether the 
Hamiltonian might be substituted by an integrable version without affecting the exponents. Examples 
are known where such a trick would not work: for instance, the higher odd spin integrable SU (2) spin 
chains flow to level k = 2s WZW models, while the chains with Heisenberg type couplings flow to 
level fc = 1. But examples arc also known where the trick does work. Indeed, with unitary groups 
and finite dimensional representations, group symmetry on the lattice plus criticality implies that 
the continuum limit is a WZW model, the stablest of all being k = 1 [18]. Therefore, a spin chain 
with SU{n) symmetry, if critical, will generically be in the same universality class as the integrable 
Sutherland chain, that is the level one SU (n) WZW model. Criticality is however not guaranteed in 
general. For instance, if one takes SU{n) and the alternance of n and n representations, the chain 
whose Hamiltonian is the Casimir has a gap, while the integrable one is gapless, and in the universality 
class of the WZW model [19]. 

Our purpose in this paper is to explore the matter further by concentrating on the spin quantum 
Hall effect where many things are known exactly. In a nutshell, the question we want to explore is 
whether the integrable 3(8)3(8)3(8)3... sl{2/l) spin chain is massless and if so, in the same universality 
class as the non integrableS 3 (8 3 (8 3 . . . Heisenberg chain, which is known to describe the physics of 
the spin quantum Hall effect. 

A positive answer would raise hopes for an integrable approach to the ordinary quantum Hall effect. 
Unfortunately, we will see that the answer is a resounding no. The reason for this is the proliferation 
of possible fixed points that appear when ordinary compact symmetries are replaced by super group, 
maybe non compact, symmetries, and sheds light on what should be directions of further research. 

We will present evidence from various sources that the integrable alternating 3, 3 chain is in the uni- 
versality class of the SU (2/1) WZW model at level 1[20]. This model has also been called OSP{2/2)_2 
in the literature [21], and has appeared previously in the study of Dirac fermions in a random gauge 
potential[22]. In contrast, it is already known from [16] that the non integrable alternating 3, 3 Heisen- 
berg chain does not scale to a WZW model, but to a new kind of theory, where the currents have 
logarithmic OPEs, and the symmetry is not of Kac Moody type [23]. 

It is important to stress here that the WZW model is a very different theory from the one [16] 
associated with the Heisenberg Hamiltonian; in particular, it has a continuous spectrum of critical 
exponents, which is not bounded from below. Trying to infer generic physical properties from the 
integrable chain would be a considerable mistake in this case, and presumably in the ordinary quantum 
Hall case as well. 
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Figure 1: Graphical representation of the representation [0, 1/2]. 

The paper is organized as follows. In the next section, we gather information, obtained from the 
literature as well as our own calculations, about the WZW model. In sections three and four, we 
present the results of a numerical and analytical study of the Bethe ansatz equations of the integrable 
spin chain, comparing our results with the expectations for the WZW model. Section five contains 
elements of extension to the case of higher level. Section six goes back to the spin quantum Hall 
problem. Technical details are provided in the appendices. 

2 First considerations on the WZW model 



2.1 Algebraic generalities 

We start with a short review of the "base" sZ(2/l) algebra. It contains a sub sl{2) with generators 
J^,Jq, an extra U{1) generator Bq, and two pairs of fermionic generators Vf^ and . Its repre- 
sentation theory is complicated, although this is one of the best understood cases in the super Lie 
algebra literature [24], [25], [26], [27]. Typical representations are characterized by a pair of numbers tra- 
ditionally called b,j (and denoted [b,j] in what follows). Here, j is a sl{2) spin quantum number and 
takes values 0, 1/2, 1, . . ., whereas b can be any complex number. Typicality requires b ^ ±j. Typical 
representations are irreducible. Their dimension is 8j, and they decompose into sl{2) representations 
with spin j,j — 1/2, j — 1/2, j — 1 with, respectively, b numbers b,b + 1/2, b — 1/2, b (for j — 1/2, 
the value j — 1 is omitted). We introduce the notation to denote an sl{2) representation of spin 
j for which the other u{l) generator takes the constant value b, so the former decomposition reads 
p'j © p'jt.\^/\ ffi Pj_i/2 ® P'j-i- "^^^ Casimir is diagonal in these representations, and proportional to 

A particularly important representation in the WZW model will be the four dimensional [b — 0, j — 
1/2] which is represented graphically in Fig. 1. It is the defining representation for osp{2/2). 

There are many more representations. Simple atypical (that is, irreducible atypical) occur when 
b — ztj, and have dimension Aj + 1. Examples are the three dimensional representations [1/2, 1/2] and 
[—1/2, 1/2] that we will use to define the 'Hilbert' space for the quantum spin chain. Representations 
contain only multiplets p'j- and p'^'^J/j, so [1/2, 1/2] contains and pj. Similarly, representations 



-J, j] contain only multiplets with p""* and p^_|[^2' 



— 1/2 — ■ 

-1/2, 1/2] contains p^^2 and 



While typical 



representations have vanishing superdimension, simple atypical ones have superdimension equal to plus 



or minus one. 
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Figure 2: Graphical representation of the representation [0, —1/2, 1/2, 0]. 



The other atypical representations are not simple: they are indecomposable ^ but partly reducible, 
and form what is often called 'blocks'. The most important example is the one appearing in the tensor 
product [0, 1/2] (g) [0, 1/2] which decomposes as the adjoint [0, 1] and another eight dimensional repre- 
sentation usually denoted as [0,-1/2, 1/2,0]. This representation is the semi direct sum of [1/2, 1/2], 
[—1/2, 1/2], and of two [0,0] representations. Its su{2) content is J^i/V^ ' Po^ Po^y ^^'^ Po twice. 
A graphical representation is given in figure 2. Note that this representation has vanishing superdi- 
mension. 

It is useful to represent such an indecomposable representation by its 'quiver diagram' [28]. Denote 
by j the representation [j, j] and by ~j the representation [— j, j] (these simple atypical representations 
can be obtained by maximum supcrsymmetrization of the fundamental or its conjugate - 2j boxes in 
a Young diagram representation). Then the eight dimensional indecomposable can be represented as 



^1/2 



\ 



1/2 



(2.1) 



Here, the bottom part of the diagram means that the representation contains a singlet as an invariant 
subspace. There is then two submodules of dimension three which are not invariant but are (and then, 
isomorphic to 1/2 and —1/2) modulo the singlet, and finally, on top, a submodule of dimension one 
which is invariant modulo all the others. A generic element of the algebra, written in the basis made 
of the top, middle and bottom components of the quiver in this order, will have the form 



/ 








a 


1/2 
^1/2 





b 





-1/2 
Pl/2 




d 


e 







(2.2) 



Here, the p are 3x3 matrices, a and b are 3 x 1, d and e are 1x3, and c is 1 x 1. 

More generally, replacing by j, and 1/2, —1/2 by j — 1/2, j + 1/2, or the same with the conjugate 
of all these representations, gives the quiver diagram for indecomposables of dimension 8 (for j = 0) 

'^We may use indecomposables as short hand for non simple indecomposables, while in the mathematics literature, 
indecomposable encompasses irreducible. 
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and 16j +4 (j > 0). These representations (denoted 1/2, j + 1/2, j]) and their conjugates are the 

only indecomposables whose quiver has a 'loop'. They are projective representations (see the appendix 
for some algebraic reminders) . All the other indecomposables ^ are not projective, and associated with 
one dimensional quivers. They will be discussed soon. 

We now turn to the current algebra. First, a note about conventions. We chose normalizations 
such that the SU{2/1) WZW at level k contains a sub SU(2) current algebra at level k in the standard 
convention for the latter. Other conventions would denote this model as the OSP(2/2) level —2k model 
(the two superalgebras are isomorphic). 

The basic OPEs read: 



J+(^)J-(0) 


k 2 Q 

= ^ + - 

z z 




J3(z)j±(0) 


= ±^ 

z 






B{z)B{0) = - 


k 
'2^ 


J+{z)V-{Q) = — 
z 


J+{z)W-{0) 


_ W+ 
z 


J-(^)F+(0) = ^ 
z 


J-{z)W+{0) 


w- 

z 


j\z)V^{^) = ±^ 


J^iz)W^{0) -- 


~2z 


ypm 

B{z)V^{0) = ^ 


B{z)W'^{0) = 




W+{z)V+{0) = — 

z 


W-{z)V-{0) 


J- 
z 



y+(z)v^-(o) = -\ + ?—^ 

2^ z 

y-(^)W+(0) = + (2.3) 

z z 

and the generators of the base algebra obtained by taking the zero modes of the currents. The Sugawara 
construction then leads to vanishing central charge independently of the level k, as the adjoint has 
vanishing superdimension. The conformal weights for afhne representations [b,j] based on irreducible 
typical representations [b, j] read 

f - 

Affine highest weight representations built on simple atypical representations have conformal weight 
equal to zero. Finally, for those built on indecomposable atypical representations, Lq is not diagonal- 
izable. 



2.2 The case k = 1 and a free field representation 

The simplest example, and the one most likely to be encountered as the continuum limit of a lattice 
model, is the case k = 1. It is well known [20, 21] that this theory admits a free field representation, 

based on a pair of symplcctic fcrmions 1^1,1^2, and a pair of bosons 0, 0'. The boson (j) is compact and 
has the usual metric, the boson (j)' has a metric of the opposite sign. The propagators are 

{r]i{z,z)r]2{w,w)) = -Inlz-wp 

{(f){z)(f){w)) = — ln(z — w) 

(0'(z)0'(w)) = ln{z-w) (2.5) 

^We restrict here to cases where the Cartan subalgebra is diagonalizable, which should be the relevant one for our 
problem. 
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The central charge of the model is — 2 + 1 + 1 = 0, as required for the SU{2/1) WZW models. 
The currents admit the following free field representation: 



J+ = 




J- = 




J3 = 


-^idd) 






B = 





(2.6) 



The fermionic currents consist of two sl{2) doublets with b = ±1/2: 

W± = e75'(±'^-'^')5,72 (2.7) 

The WZW model has never (to our knowledge) been completely understood, even for fc = 1: 
questions such as the operator content or the operator algebra only admit partial answers, at best. 
Of course, it is tempting to speculate that the operator content derives somehow from the sub SU (2) 
part, which for k = 1 contains only the affinc representations at j = 0, j = 1/2. The j = 1/2 case 
corresponds, for the fundamental representation of osp{2/2) which has 6 = 0, to a conformal weight 
/i = |, and one thus expects to have, for A; = 1, at least the representations with h = and /i = g. 
But it is easy to argue that the operator content is considerably richer, and in particular exhibits 
logarithmic features. Indeed, consider the tensor product of the [0, 1/2] representation (the multiplet 
for the fundamental field of the theory) with itself; 

[0, 1/2] (gj [0, 1/2] = [0, 1] + [0, 1/2,-1/2, 0] (2.8) 

The right hand side does not contain the true singlet [0, 0], which means, jumping from tensor product 
in the base algebra to fusion product in the conformal field theory [21], that the identity field must 
have logarithmic partners, and fields of dimension h = he organized in several representations of 

5^(2/1). 

To gain further insight on this question, we can study the OPE associated with (2.8) by using the 
free field representation. The doublet in the [0, 1/2] representation can be represented by [21] 

v±=e'^^IJ, (2.9) 

where /x is the A = ^ twist field in the symplectic fermion theory, with dimension —1/8, which, added 
to +1/4, gives the desired 1/8. The two singlets meanwhile can be represented as [21] 

w±=e'^^iy'^ (2.10) 

where are the A = — 1/2 and A = 3/2 twist fields in the symplectic fermion theory. Their dimension 
is 3/8, which, added to —1/4, gives the desired 1/8. The fields in the adjoint representation are the cur- 
rents, whose expression was given earlier. Finally, for the fields in the indecomposable representation, 
they are given by 

111112 

1 (2.11) 

where the fields are arranged exactly as in the quiver for [0,-1/2, 1/2, 0] given above. All these fields 
have vanishing conformal dimension. The superdimension of [0, —1/2, 1/2,0] vanishes. On the other 
hand, it is known (see below) that the Witten index of the theory is equal to one. This means that in 
the level fc = 1 WZW model there must be at least the fields in (2.11) and a 'true' identity field, which 



e=^75e ^r/i, e'''^'''' rjidrji 
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is not embedded in a bigger indecomposable representation. The free field representations is thus not 
complete as is, since two fields '1' must be introduced. 

In fact, it is possible to build (many) other indecomposable representations associated with fields 
of vanishing conformal weight. The simplest are obtained as representations containing either the field 
r/i or the field r/2- For instance, the fields 

^2, e75(±*+'^'), e*^^'9ryi (2.12) 

form a representation one can represent by the diagram 2^ • . Similarly, the fields 

^1, e7!(±'^-'^'), e-^^'^'drj2 (2.13) 

can be represented by the diagram •'^ <— •. The module at the end of the arrow is invariant, and 
quotienting by this module amounts to erasing the dot, that is the remaining module is then invariant. 
For more general such diagrams, the dots will represent atypical representations of the type [±i,i] 
discussed earlier. 

These structures do not come alone, but are always embeddable in larger structures (in mathemat- 
ical terms, they arc not 'projective' representations). This is easily seen in the free field representation. 

For instance, introducing the field e^^^'^~'^'^ri2 and acting with the generators gives rise to • • <— •. 
where the [—1,1] submodule is constituted by the fields 

giV2(0-0')^2, e75W-<^')^2m, e-^^^'r;2, eVf (^+^')r;2,yi, e-'v^(^+^')r?2 (2.14) 

In turn, these fields can be connected to bigger representations, and one finds that the conformal fields 
of vanishing weight built 'on top' of a single fermion r?i or 772 can be arranged into the semi-infinite 
quiver 

-p -V+\ -p+l -1 -1/2 /r,ir\ 

• <— • • -i— • (2-15) 

or its mirror image. The highest weight (that is, the field annihilated by the zero modes of all raising 

operators) of — n/2 is the field e^^'^~'^ ^2 for n even, e^^*^""^ ^ for n odd, and similarly for +n/2. 

Of course, since we are at level k = 1, these higest weights for n > 1 are not aSine highest weights, 
that is, their OPE with the SU{2) currents gives poles of degree higher than one. For instance for 
n = 2, 

j-(^)e^V2(0-0')^2(w) = , ^ : e-'^^(^)e'^(^-^')r72(«;) : (2.16) 

(z — w)'' 

Therefore, applying modes with positive indices will produce fields with negative dimension. In this 
example we get the field e"*^"^ 772 of dimension —1. This field belongs to the representation [b = 
—3/2, J = 1/2], for which e^^^^'^~^'^'^dr]2 is an afiine highest weight. The pattern generalizes, and one 
finds that for any n > 1, the fields at weight /i = in the semi-infinite quiver for node —n/2 are affine 
descendents of the field 

dT^2.-.d-^-^m (2.17) 



exp 



-J=(0-(2n-l)0') 



which is the highest weight of [— n -|- 1/2, 1/2], with h = —n{n — l)/2. Concatenating [—1/2, 1/2] and 
[0,0] in a four dimensional representation we will still denote (abusively) by [—1/2,1/2], the set of 
conformal fields associated with the semi-infinite quiver seems to be made of all afiine highest weight 
representations with base [b = —n + 1/2, j = 1/2], n = 1, ... ,00. Similarly for the mirror image, we 
get all the representations [b = n— 1/2, 1/2], n = 1, . . . , 00. 

Note that as we increase the value of n we get larger and larger negative dimensions, and the 
semi-infinite quiver corresponds to a theory whose spectrum is not bounded from below. This feature 
is rather characteristic of supergroup WZW models. It does not seem to make physical applications 
impossible however, as we will discuss below. 

The algebraic structure of the representations of the affine algebra associated with the semi-infinite 
quivers is very interesting of course, and deserves much further study, beyond the scope of this paper. 
We shall refer to them, whenever necessary, as 'infinite blocks' representations of the current algebra. 
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Note that in a physical theory, left and right sectors have to be combined. While we have treated 
the fermions so far as chiral objects, they are not, and when completing with the left moving parts of 
the boson operators, the left moving part of the fermions should also be included. 

Let us now backtrack a bit and wonder what of this algebraic zoo will appear in the WZW model. 
First, observe that the free field representation of the currents uses only the derivatives of the fermions, 
which arc well defined conformal (chiral) fields. One might have hoped that maybe a theory can be de- 
fined that does not use the fermion fields themselves (and thus is based only on the "small algebra" [29]). 
However, consideration of the four point functions and the Knizhnik Zamolodchikov equation seem to 
exclude this possibility, as the identity clearly must have logarithmic partners, and the fermions them- 
selves are necessary to reproduce this feature in the free field representations [21]. One could then 
have hoped that fields of dimension zero are reduced to the true singlet and the eight dimensional 
indecomposable [0, 1/2. —1/2, 0]. This seems also incorrect. For instance, it is known that the fusion of 
the twist fields with dimension h = —1/8 and /i = 3/8 in the symplectic fermions theory does produce 
the fermionic fields themselves [29]. Also, we shall see later that modular transformations of characters 
suggests the presence of more fields of weight zero. If so, then the non projective indecomposables 
must appear, and then there is no reason to restrict them. We in fact speculate that the infinite 
blocks defined previously appear in the WZW model at level one. This in turn implies apppearance of 
arbitrarily large negative conformal weights. 

Indecomposable representations of s^(2/l) lead to affine highest weight representations with van- 
ishing conformal weights. To have a non vanishing conformal weight for an affine highest weight 
representation requires dealing with a typical, irreducible representation of sl{2/l). There, the pres- 
ence of a sub SU{2) at level one presumably allows only spins j = 0, 1/2 to occur. So a minimal guess 
for the WZW model would be to add to the infinite blocks just discussed, the affine highest weight 
representation based on [0, j = 1/2]. But we have already seen appearance of the representations 
[— n-|- 1/2, 1/2] in relation with the infinite blocks. Thus it seems to us that the spectrum of the WZW 
model should at least include 

oo 

[M]®[M] [-n + 1/2, 1/2] ® [n -I7271/2] (2.18) 

n= — oo 

2.3 Relevance to lattice models? 

The appearance of arbitrarily large negative conformal weights is related with the fact that the invariant 
metric on the supergroup is not definite positive, and thus the functional integral ill defined - a feature 
directly at the origin of the 'wrong sign' in the propagator of the free field (p' . This might lead one to 
think that the conformal field theory we have tried to define through the free field representation of 
the SU(2/l) current algebra cannot in any case appear as the continuum limit of a lattice model, for 
which (sec below) one can easily establish the existence of a unique ground state;. But things are more 
subtle. It might be for instance that the arbitrarily large negative conformal weights correspond 
to non normalizable states, and thus are absent from the spectrum and the partition function. This 
would be the situation for a lattice model whose continuum limit is a non compact boson [31] . It might 
also be that the arbitrarily large negative conformal weights do correspond to normalizable states, but 
that the correspondence between the spectrum and the partition function is more complicated than 
for unitary conformal field theories say, and proceeds through some sort of analytic continuation. This 
seems to be the case for the (3^ system for instance [32], where instead of the naive characters, the 
modular invariant partition function is made of character functions, which encode the arbitrarily large 
negative conformal weights in a subtle way. We will argue that this is the situation in the SU{2/1) 
model as well. In other words, the spectrum deduced from the lattice model may be related to the 
spectrum of the conformal field theory only through some sort of analytical continuation. We will 
discuss below a very naive version of this continuation. 

Finally, and getting a bit ahead of ourselves, one might recall that Zirnbauer (sec [33] for a thorough 
discussion) has suggested to replace the target space from (in this case) S'L''(2/1) with its indefinite 

^Note that these axe not necessarily excluded on "physical grounds'. It might well be that in some systems, these 
axbirarily large negative conformal weights do describe physical observables, like powers of the wave function [30] . 
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metric by a non compact target space with a positive metric - a Riemannian symmetric superspace. 
The resulting functional integral would then be well defined, and the spectrum bounded from below. 
It is not clear to us how one would study this model in practice. But it is tempting to conjecture that 
its spectrum would coincide with the analytical continuation of the naive spectrum deduced from the 
algebraic study of the current algebra, that is the spectrum observed in the lattice model. If so, the 
question of 'what is the continuum limit of the lattice model' may have, in a way, two answers. This 
obviously requires more thinking. 

To proceed, we now turn to characters and their modular transformations. This will give us a handle 
on the mysterious analytical continuation at play, and allow us to make contact with the conjecture 
(2.18). 

3 Some results about the WZW model from characters 

3.1 Integrable characters. 
3.1.1 Ramond sector 

At level A; = 1, four admissible characters [34] have been identified. In the Ramond sector (periodic 
fermions on the cylinder) one has the field ^ = 1/2,5 = with character 

R,I ( \ „l/8,l/2 TTii?/^ ^ ,A \^ „2a^ +0/^20 ^-2a-l-i (n i\ 

Here we use the definition ^ = Tr qLo-c/24^2nraJ^ ^2%TnyBo ^ j3 ^j^j j^^ being the zero modes of the 
field (the third component of the spin) and the B field. Normalizations are such that for the highest 
weight state in the /i = 1/8 R sector, = j = 1/2, Bo = b = 0. We have set z = e^''^" and C = e^^""" . 
Finally, 

_ (1 + z''\y^q-){l + z-V\V^g"-i)(l + zV^r ^/^g")(l + z-'/^-^^q-^) 

^ (9'^'0-ll (l-g")2(l-zg")(l-z-V-l) ^^■'^> 

It is convenient to rewrite this character as [34] 

Xhdi/sir, i^) = ^ [xo{t, (t)xi/2{t, u) + Xi/2(r, (j)xo{t, u)] (3.3) 

where the x's in (3.3) are characters of the SU {2)i WZW model, and are given by the following well 
known expressions [35] 

XiMr,a) = ^^Y.^^'-^'^'^z'^^y^ (3.4) 



(note that these characters are even functions of cr, as expected by symmetry of the Jg spectrum), and 
we have set z = e^^^^; similarly we set ( = e^"^". The specialized Ramond character a = v = for the 
h = l/8 field reads therefore 

Xhdi/si^) = ^XoXiMr) = 4gi/« (1 + 5g + 22q^ + . . .) (3.5) 

where we recognize the overall multiplicity 4, coinciding with the dimension of the representation 
[0, 1/2]. The latter decomposes as one spin 1/2 and two spin SU{2) representations. 
The Ramond character for the identity field is considerably more complicated [34] : 

X^io" = F-ir, a, .) ^ ^-'^ ^ ^ ^,-y2^;!;^^l ^,-^2^/2^ (3-6) 
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The specialized character turns out to be finite, despite the divergence of F^: 



R.IV 

Xh=o 



1 „a °° „2a^+o 

0=1 a=l ^ ' 



1 + 8g + 24g2 + . . . (3.7) 



We thus sec that in the Ramond sector, all conformal weights arc of the form h ^ Mod integers or 
/i = I Mod integers. Setting a ~ 1 gives us Tr(— 1)^ (a supertrace) and evaluates what is called the 
super character; we check that x'h^i/si''') ^ while Xh^o^ {t) = 1- 



It is possible to write the specialized character in a more compact form using the function [36] 

1 



One finds then 



2 2 



« = ^2^J (^) + 2X0Xi/2/i3(r) (3.9) 

AO "T Ai/2 

3.1.2 Neveu-Schwarz sector 

In the NS sector, the fermions have antipcriodic boundary conditions, and the SU{2/\) symmetry is 
broken down to a sub SU{2) x U{1). The lowest energy state in this sector has negative conformal 
dimension, leading to an effective central charge equal to Cg// = 3. The corresponding character is 



Xh=-i/8i'^'^'^) = [xi/2(T,cr)xi/2(T, i^) + Xoir, a)xo{T, u)] 
The specialized characters reads 



(3.10) 



NS. 



-V8W=9-'/'(l + V/^ + ...) (3.11) 



so half integer gaps appear in this sector. 

The other character is again more complicated 



1 — zq 

aez 

where. 



,2a+l 



Xh=i/4l^'«^''^J g 2 ^ Kr,(J,v) yq z ^ ^„+i/2^i/2^_i/2)(i ^ ^a+1/2^1/2^1/2) 

(3.12) 



_ ^ (1 + ^V2^V2^"-V2)(l + ^-l/2^1/2gn-l/2)(^ ^ ^l/2^-l/2gn-l/2-)(i ^ ^-l/2^-l/2^n-l/2-) 

(l-g»)2(l-zg")(l-z-ig»-i) 

(3.13) 

The specialized character is again finite, despite the divergence of F'^^: 

XLY/IW = hz{T) [xo(t)2 ^ xi/2{rf] (3.14) 
and like Xh=-i/8' exhibits half integer gaps. 
3.2 The remaining operator content 

It is crucial to realize that the characters we wrote down all have complicated or unusual modular 

transformations. A possible strategy to proceed is then to try to supplement these characters with 
others in order to obtain a finite dimensional representation of the modular group, and maybe modular 
invariants. 
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3.2.1 A continuum above /i = 1/8 in the R sector 



• Transformation of x^^-i/s 

Consider for instance the Neveu Schwarz characters. They correspond to fermions being antiperiodic 
in both directions. Aternatively, this corresponds to an antiperiodic chain along the space direction, 
while one takes a trace along the (imaginary) time direction. Imagine instead taking a supertrace along 
the imaginary time direction, or setting cr = 1 in the expression of the characters, giving rise again to 
a supercharacter. The specialized supercharacter for /i = —1/8 reads 

2 2 

sNS,I / \ '^O ~ Xi/2 , , /Q 1 r\ 

Xfc=_'i/8W = (r) (3.15) 

Recall now the modular transformation of the basic 5C/(2)i WZW characters 

) i-l/r) = f ? )ir) (3.16) 

Xl/2 J \ V2 72 J \ ^1/2 / 

Meanwhile the i] fimction obeys ?](— 1/t) = \/ —iTTjir). So the only way to write the modular transform 
of this SU {2/1) (super) character as a sum of powers of q is to introduce an integral representation of 
the T] factor 

1 



and thus 



oo 



rfae"'™ (3.17) 



J''-^=^hd,/sJ daq"/' (3.18) 

Now observe that under modular transformation the boundary conditions which determine the super 
NS character turn into boundary conditions where the fermions are periodic in the space direction, and 
one takes a trace in the imaginary time direction. Hence the right hand side of (3.18) should appear in 
the generating function of the spectrum in the sZ(2/l) symmetric sector of our chain if the associated 
conformal field theory contains the Neveu Schwarz characters discussed previously. This means there 
should be a continuous component stcu-ting at the h = 1/8 value. 

To find out about this continuous component, we can keep track of the Bq and Jq numbers by 
keeping the parameters tJ, u: 

xll%{T, <T, u) ^ Tr [{-Ifq^^'^z^k''''] = ^ [xo(r, a)xo{r, u) - Xi/2{r, a)xi/2(r, i^)] (3.19) 

The modular transformation of general SU (2) characters reads 

^° ) (-l/r,e/r) = e^-^V^^ f ^ ] ( ) (r,9) (3.20) 

Xl/2 J V 72 71 / V Xl/2 J 

and thus 

X^=!'V8(- C7/T, v/t) = 1 iA^^+^^)/2r [^^(^^ a)xi/2(T, u) + Xi/2(r, a)xo{r, i^)] (3.21) 

Due to the indefinite metric for the SU {2/1) supergroup, the overall factor in modular transformations 
should however not be e*'^^'^ )/'^'^^ but rather e™^" ~" ^Z^'" (since the Casimir is proportional to 
j'^ — 6^). We are then left with an overall factor e^'"" 1'^ whose interpretation is slightly delicate. What 
we propose is to write 

V-iT V-1 J-oo 
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The phase factor is ambiguous and depends on a choise of cut. As for the integral, is is divergent in 
the physical situation where Im r > 0. If we assume however that this integral is defined by some sort 
of analytical continuation from the case Im t < 0, wc sec that wc can formally interpret the leftover 
exponential in our modular transform (3.21) as resulting from an integral over fields with conformal 
dimension and Bq charges: 

/i = -^6^, Bo = 6, -oo<6<oo (3.23) 
Moreover, the right hand side of the modular transform is, for each value of 6, of the form 

-^^^-'^^"^^ [xo(r,a)xi/2(r,i^) +Xi/2(r,a)xo(r,z^)] (3.24) 

which is exactly the character for the affine representation with base [6, 1/2] (integrable character in 

class I in the language of [34]). It thus seems that the model should contain in fact, not only the 
representation [0, 1/2], but the continuum [6, 1/2] with h real. This is at odds with the first section 
where we speculated that [— n + 1/2, 1/2] only apppeared. See below for a discussion of this issue. 

While the conformal weights associated with those representations become arbitrarily large and 
negative, the analysis of the characters and modular transforms suggests that the finite size spectrum 
is given by a sort of analytical continuation, whose result, by the mechanism discussed above, is to make 
the continuum appear above the 1/8, not below. Note also that while the Bq numbers in the chiral and 
antichiral sector appear continuous, the lattice analysis anyway only identifies the combination L + R, 
which is made of integers. 

• Treinsformation of XhLifi 

We can perform the same analysis for the specialized supercharacter for /i = 1/4: 

sNS.IV / \ o XOXl/2 , , / 2 2 \ /n oc\ 

Xh=i/i W = -2-2— 2- + /j3(Xo-Xi/2) (3-25) 

XO "1" Xi/2 ^ ' 



The function /is obeys 



/.3(-l/r) = -Mr) + i£c/a^£;i^ (3.26) 



Thus 



2 2 

sNSJVi -, / N Xl/2~Xo, s , „ / X 

Xft=i/4 (-V^) = J .o ij) + 2xoXi/2(t) 

Xl/2 Xo 



1 g«'/2 
-hair) + — / da- 



r]{T) J_^ 2cosh7ra 



(3.27) 



where we again obtain a continuous spectrum starting at /i = |. This time however, instead of a flat 
measure over all the reals, one has a measure proportional to l/cosh7ra. Note we can write 

8NS,IV^ 1 / N RJV^ ^ , R,I , Q ' /o oo\ 

Xh=i/4 (-V^) = -Xh=o (^)+XhU/8 J ^" 2cosh7ra ^ ^ ^ 

To find out about the spectrum of Bq and Jq values is considerably more complicated, and the 
tired reader may want to jump to formula (3.44). The following calculations rely on the paper [37]; 
one starts by writing 



Xfj}7{r, a, V) = _e^-/2e^-F«^^(r, a, u) 
where 



T + a V -T-a V 

i^4(r,^,--)-i^4(r,^— ,--) 



(3.29) 



_ -pr (1 - zl/2^V2^"-l/2)(;^ _ ^-l/2^1/2^n-l/2)(i _ ^l/2^-l/2gn^l/2) _ ^-l/2^-l/2^n-l/2) 

~ (1 - g")2(l - ^;5")(1 - 0-lg"-l) 

(3.30) 
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It is useful to rewrite this prefactor as 



e,,{T,a)q-^/W{q) ^ ' ' 



giTrm lr-\-2i7rmla 



The K are Apell functions 

KiiT,a,.)=Y. (3-32) 
and the 6 are usual theta functions 

oioiq, ^) = n (1 + ^"'«'") n (1 + ^9") n (1 - 9™) 

m>0 m>l m>l 

^11 (9. = n (1 - '^"'9") n (1 - ^9") n (1 - (3.33) 

m>0 m>l m>l 

Modular transformations of the theta functions are as follows: 

^io(-l/r, /x/r) = x/=i^e"^^^eio(r, ^ + 1/2 - r/2) (3.34) 
As for the Appell functions it is much more complicated. One finds [37] 

Ki{-1/T,a/T,iy/T) = e"'^''"-''"'>/^Ki{T,a,iy) 
l-i 

+T e''^'('"+«^/')'/^$(ZT, lu - aT)eoo{lT, la + ar) (3.35) 

a=0 

where 6oo{q,z) = Oi(,{q,zq~^/'^) and 

'^(r,,) = --^ - \ I ,,e-- -^"^(--^(l±_^;/-)) (3.36) 
2\/—iT i J smh(7ra;v— «t) 

Putting everything together gives 

X'^=f/f (-1/t,(7/t, vjr) = ie--(-'-')/2-e2Wrg-Wr 
OiojT, -(g + t/)/2)^io(T, -(g - u)/2) 

[K4{-1/t, (<7 - l)/2r, -Z//2T) - K4{-1/t, {-a - 1)/2t, -!//2r)] (3.37) 
Using the general Appell transform formulas gives 

K^i-l/r, {a - 1)/2t, -iy/2T) = Te'''^'''-''^^/^e'''''e-^''"'/^K4T, {a - l)/2r, -i^/2) 
i-i 

+T e*'^('^-i+"^/=^)'/^$(4T, -2u - ar)6»(4r, -2 + (2<t + a)r) (3.38) 

0=0 

We recognize the expected global factor e^'^(<^^-'^^)/^'r _ Also, observe that 

0io{q,xy)6io{q,xy-'^) ^ eiojq, x-^y''^)eia(q, x-'^y) 

eii{q,x) Oiiiq,x-^) ^ ■ ' 

while 

pR. ^ . _ Oiojr, {a + ^)/2)0io(r(a - u)/2) 
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We also observe that K4{t, v, fx) = K4{t, v ± 1, At)- Therefore, the contribution to (3.37) coming from 
the Ki term in the transform (3.38) gives 

First term = -e^^'^'-'^'^/^r^W^^^ ^^.Al) 
The second term reads, after making the 6 function exphcit, 

3 

-F«(t, cr, i/)ie*''('^'+''')/2^ Y^{-IY^{At, -2v - aT)q''' ^^z"/^ 9^'"'+""' (z^"" - z-"^""-") (3.42) 

a=l m 

We recognize the prefactor of the $ function as the formal extension of the Ramond character formula, 
initially valid only for j = 1/2, b = 0,h = 1/8 to j = a,b = 0,h = Xh=a^/s- However, setting 

Xa = g"'/^^"/^ J2 9^"'+""" (z^" - 0-2™-") (3.43) 

m 

it is easy to show that Xa = Xa-4 = —X-a- Therefore, Xq = X±2 = 0, and Xi = — X_3. Taking 
advantage of these symmetries leads us to a considerably simplified expression 

+^-dysir,^,.)e-i^'^-^'y^^ [ cosh(2...v^/r) 
' J 2 cosh(7rxv— «t) 

As a check we can set z/ = in this expression. Using the invariance of hs{T) + /i3(— l/r) we see that 

r?(r) J 2cosh7ra ^ ^ ^ ^ 

is invariant, and thus 

/ dx ^= = / da— ^ (3.46) 

J 2cosh7rxV-»r j 2cosh7ra 



which alows us to recover expression (3.28). 

We now go back to the more general expression. Crucial is the fact, as mentioned earlier, that the 
exponential has the v"^ term with positive sign instead of a negative one. We thus write in the integral 

X = y\J^ and thus formally obtain 

V T J coshiTry 

dy (3.47) 

cos Try 

Using the formal integral representation of the exponential we write this as 



COS Try 



/■ 
dbdye-''^''^^ ^ 

= / d6e2*'^''''e-"^'''^— (3.48) 

J cos TTO 

Finally, we thus rewrite (3.44) as 

xSr(-l/r,a/r,./r) = -e-(^^-^^)/2^x^io^(r, a, .) 
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This corresponds again to representations [b, 1/2] but this time with a measure that is not flat, but 
proportional to 1/ cos nb. The interpretation is discussed below. 



• Transformation of Xhl-i/s ^h^i/i 

Finally, wc can also consider the ordinary NS character. Wc know that the boundary conditions 
corresponding to it character are invariant under the modular transformation. But using the expression 

2 2 

NS,I , X ^0 + Xi/2 , , , . 

XhJ-i/si.-r) = {r) (3.50) 

we see that the modular transform has a continuous component starting at the ground state of the 
sector, namely /i = — 1/8: 

xL":'i/8(- Vr) = ^xL":fi/8 W = {xl + Xl/,) (r) / da^-^ (3.51) 

Similarly 

/oo „a^/2 



A final note: in the supersymmetric literature (the model has some relations \.o N = A. theories 
[36]), the subset with discrete quantum numbers is refered to as massless, the one with continuous 
quantum numbers as massive. 



3.3 Interpretation 

Comparing the free field analysis of the WZW model with the analysis of the modular transformations of 
basic characters suggests that the spectrum of arbitrarily large negative dimensions admits 'regularized' 
expressions for characters, and thus partition function and spectra of lattice Hamiltonians. For these, 
there is a single ground state, and the negative dimensions are 'folded' back into a continuous spectrum 
of positive dimensions. 

The slightly surprising thing is that wc find in this way indications of a continuous spectrum 
based on [6, 1/2] while we initially expected a discrete one [— n+ 1/2, 1/2]. There are stronger reasons 
for a discrete spectrum of b than the ones given in the first sections. For instance the free field 
representation of the generators selects a particular radius for the non compact boson, and vertex 
operators corresponding to charges continuous charges b will in general be non local with respect to 
these currents! Of course it might be in fact that the free field representation we have used is good 
only for a subset of fields, and docs not allow one to explore the full operator content. But it could 
also be that the continuous spectrum we read from the modular transformations is an artifact from 
the (rather uncontrolled) regularization at work. A hint in this direction comes from the integral 

J 2 cosh 7ra 

which, after continuation to the imaginary axis, has poles precisely at the values b = —n + |. It seems 
diflacult to conclude on this issue without more work on the conformal field theory side. 



4 Analysis of the lattice model 

The next step is to extract as much information about the lattice model as possible and compare it 
with our proposal for the continuum limit. Although we have mostly talked about a quantum spin 
chain so far, integrability allows one to consider it as part of a larger family of commuting transfer 
matrices related to an integrable vertex model. 
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4.1 Definition of the integrable vertex model 



The vertex model is constructed from ii-matrices acting on tensor products of spaces carrying the 
fundamental representation 3 of the graded Lie algebra sZ(2/l) and its dual 3. The i?-matrices satisfying 
Yang-Baxter equations (YBEs) have been constructed in [38, 39, 40] (see also [41]) and read 



V 

1 + -0 

V 



i?33(u) = 1--P 
V 



(4.1) 



Here V and O are the graded permutation and monoid operators of s/(2/l), respectively (note that 
the action of the monoid O on 3 (S) 3 differs from that on 3 (S) 3). As a consequence oiO^ =0 we have 

it!33(-«-l)i?33(«-l) = l. 

Prom these ii-matrices we construct two families of monodromy matrices (labelled by the spectral 
parameter ?;). The monodromy matrices act on the tensor product of a "quantum space" (3 $5 3)^^ 
(the Hilbert space of the superspin chains considered below), and a three-dimensional "matrix space" 
carrying a representation 3 (3) 



1 



X)R^i^iv)Rf^iv 



1 



X)---Rg\v)R';^\v-l-X) , 



i?«(. 



1 + A)i?,^y (i;)i?|^ (?; - 1 + A)i?,|^ (v) ■ ■ ■ i^^g^ {v - 1 + A)i?_^f ^ (v) 



(4.2) 



Here the parameter A can be chosen freely without affecting the essential properties (symmetries and 
locality of the interactions) of the models considered in the following. Variation of A allows one to tune 
the coupling between the two sublattices of 3 and 3 spins, respectively. Later we will focus on the case 
of A = which corresponds to the most symmetric coupling. As a consequence of the YBE 



the corresponding transfer matrices 



i?33(u - V) niu) -1-X)=T-3{V-1-X) n{u) R^-^iu - V) 



Ta{v) = Stro (Taiv)) , 



3,3, 



(4.3) 



(4.4) 



commute with each other for arbitrary values of the spectral parameter v. Here stro denotes the 
supertrace over the matrix space. 

By staggering the transfer matrices (4.4) one can construct the integrable vertex model shown in 
Fig. 3. Arrows pointing upwards ot to the right (downwards or to the left) indicate that the corre- 
sponding link carries the fundamental 3 (the dual 3) representation of sl{2/l). The vertex weights 

3 3 



3 3 



Figure 3: The vertex model. 

corresponding to the various arrow configurations are shown in Fig. 4. The transfer matrix of this 
model includes two layers in the vertical direction and is equal to t{v) = T3(v)t^{v). 

The vertex model constructed in this way was suggested as a description of the spin quantum Hall 
effect in Ref. [40]. The notations of [40] are obtained by setting v = —2x and A = 2u — 1, where x and 
u are quantities defined in [40] . 
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[^33(^-1)]^^ 



Figure 4: The vertex weights. Arrows pointing upwards ot to the right (downwards or to the left) 
indicate that the corresponding hnk carries the fundamental 3 (the dual 3) representation of si (2/1). 



The Hamiltonian of an integrable supersymmetric spin chain with local (nearest and next-nearest 
neighbour) interactions only is obtained by taking the logarithmic derivative of the two-layer transfer 
matrix at w = 0'' 



(4.5) 

where a{v) = 1 — 2/v ^. The Hamiltonian (4.5) is not hermitian and its eigenvalues are in general 
complex. However, for small systems we find that the ground state and low-lying excited states have 
real eigenvalues. We assume that this continues to hold true in the thermodynamic limit. Under this 
assumption we find that in the thermodynamic limit low lying excited states can be understood in 
terms of a conformal field theory and we can extract the central charge and scaling dimensions from 
the finite-size energies for the ground state and low-lying excited states. 

The Hamiltonian (4.5) contains a piece proportional to the sl{2/l) invariant product [42] Sj.Sj+i 

where Sj are sl{2/l) generators on site j, and the contraction defines the invariant form. In addition 
to this Heisenberg coupling however, our Hamiltonian (4.5) contains more complicated terms involving 
three neighbouring sites. 

The eigenvalues and eigenvectors of the transfer matrices were constructed by Links and Foerster 
in [38] by means of he graded quantum inverse scattering method [43, 44] . As usual different choices 
of reference state lead to different forms for the Bethe ansatz equations [45, 46, 44, 47, 48, 49, 50, 51, 
52, 53]. In the following we will analyze the Bethe ansatz equations for the choice of grading "[010]" 
which read 

!^Y = Yl^^l^Jl±l, j = l,...,N 
l^ll^)' ^TTln^lJlll^ a=l,...,M. (4.6) 

Equations (4.6) coincide in fact with the BAE describing the spectrum of the so-called "Quantum 

Transfer Matrix" of the ID supersymmetric t J model see Appendix D. 

Each solution of (4.6) parametrizes an eigenstate with quantum numbers B — (N — M)/2 and 
= L - (M + A^)/2 (we use the notation B = Bq + B° and = + J^). The eigenvalues of the 

"^A similar Hamiltonian was considered in [38] from the point of view of introducing impurities into a t-J model. 
^As usual the overall scale of the Hamiltonian is at our disposal. 
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transfer matrices (4.4) are [38] 



TV M 

A{v) = [a(u)]'^ Yl a{l + iuk - v) + [a(l + A - f)]^ [] a{v - ij/s) 

fe=l /3=1 

TV M 

- a(l + iuk - ^7/3) ) 

k=l /3=1 

(4.7) 

M TV 

A(?;) = [a{v)]^ Yl a{iji3 + l-X-v) + [a(l - A - v)]^ [| a(i; + A - m^) 

/3=1 fe=l 

TW TV 

- J]^ a{iji3 + 1 - X - v) Yl_a{v + X - iuk) ■ 

/3=1 fc=l 

The eigenvalues of the Hamiltonian (4.5) are obtained by taking the appropriate logarithmic derivatives 
of the eigenvalues of the transfer matrix 

Similarly the eigenvalues of the momentum operator are obtained are 

In the definition of the momentum we have taken into account that the underlying translational sym- 
metry is by two sites of the lattice. 

4.2 The ground state of the lattice model 

In order to gain some insight in the structure of the spectrum of the Hamiltonian (4.5) we have 
diagonalized the transfer matrix and the Hamiltonian explicitly on small lattices. For L = 1 and 2 (i.e. 
2 and 4 sites) this was done analytically and for X = 3, 4, 5 (corresponding to 6,8,10 sites respectively) 
numerically. The results of this analysis arc summarized in Appendix B.2. 

Perhaps the most striking result of this analysis is that the space of true singlets (i.e. states 
annihilated by all the generators of the algebra, which are not the image of any other state under the 
action of one of these generators) appears to be one dimensional. This is not entirely obvious: although 
the superdimension of the tensor product is always equal to one, it is a priori conceivable that there 
are more singlets, and that the superdimension one results from cancellations. But the result can be 
established rigorously (see Appendix A). Note however that one can always combine an invariant (i.e. 
annihilated by all generators) state that is not the image of any other state with another invariant 
state, and still get a state that is not the image of any other state. In that respect, although the 
singlet appears only once in the decomposition of the spectrum into representations of sZ(2/l), the 
exact formula for this singlet depends on the Hamiltonian. 

Based on our results for small systems (see Appendix B.2) we conjecture the following form for the 
eigenvalues of the transfer matrices ts{v), 73 (w) of this singlet state 

^ \ v{v-l-X) J ^ \ v{v-l + X) J ^ ^ 

Interestingly, there exists no non-degenerate solution of the BAE (4.6) giving rise to these eigenvalues. 
Applying a twist in the boundary conditions for small systems and studying the evolution of BA roots 
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as this twist goes to zero we were able to verify, however, that the singlet is in fact described by the 
degenerate solution Uk = A, 7^ = 0,k,a = 1, . . . ,L which clearly gives Ai{v) from the general form 
(4.7). This finding coincides with the observation in [54] regarding the distribution of BA-roots for 
the ground state of the quantum transfer matrix for the t-J model in the thermodynamic limit. The 
eigenvalue of the full transfer matrix corresponding to the singlet state is 

A{v) = Miv)Miv) = .>-l-A)(.-l + A) ) ■ ^'-''^ 

For A = this state is the ground state of the Hamiltonian with energy 

Eo = -AL . (4.12) 

The energy of this singlet is proportional to L without any finite size corrections which implies that 
the central charge of the conformal field theory arising in the scaling limit of the lattice model is zero 

c = . (4.13) 

We now turn to a discussion of Bethe ansatz results on the low-lying excitations of the lattice 
system. We present a rather short discussion here, and refer to Appendix B for the details. 



4.3 Excitations 

Based on an analysis of the analytical properties of the bare scattering phase shifts appearing on 

the right hand sides of the Bethe ansatz equations (4.6) we are able to classify the possible root 
configurations, so called strings, arising in their solution (see Appendix B.l). We identify the root 
configurations giving rise to low energy excited states by considering small systems. We find that 
many low-lying states are built from two types of complexes involving one spectral parameter u and 
one spectral parameter 7 each (the n = 1 "strange string" configurations, see (B.3) and (B.4)) 

c-i i (-■) i * 

type - : ' = X - - + e- , T = x + - + e_ , 

\ \ ^4.14) 

type + : vS'^^ = a; -I- - -I- e+ , 7^"'"' = x — - -I- . 

A particular class of low-lying excited states is obtained by combining A'"_|_ strange strings of type 
+ with A''_ strange strings of type — 

"1"^ = - J + e^,- = ' i = 1, • • • , , 

i * (4-15) 

= + I + = (71+^)* , fc = 1, . . . ,iV+ . 

Our numerical results indicate that the corrections e± to the "ideal" string solutions are of order 
1/L. Each of the two complexes (4.14) constitutes an acceptable, non self- conjugate, solution of the 
Bethe ansatz equations (4.6). Such solutions have been encountered in other studies of non unitary 
Hamiltonians, like the s/(3) integrable spin chain [55]. 

To make contact to the analysis of the WZW model we have to determine the spectrum of scaling 
dimensions for the field theory describing the scaling limit of the lattice model. Conformal invariance 
implies that the dimensions can be extracted from a finite size scaling analysis of the low-lying energies 
of the latter. Before evaluating these energies in large, finite systems it is very useful to consider an 
infinite volume first. In the thermodynamic limit we can neglect the deviations from the ideal string 
solutions as they do not contribute to the leading 0(1) behaviour of the energies. Then, using (4.15) 
in (4.6) and multiply the equations for the components of a strange string we arrive at the following 
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set of equations involving only the real centres of the strange strings 




y3-Xk + l\ T-r / Vj -yk + 2t 



2 



(4.16) 



fc=l ' fc=l 



The algebraic equations (4.16) can be turned into coupled integral equations for root densities in the 
limit L 00, N± oo, ^ = n± fixed (see Appendix B.3) 

p+{x) + p'l{x) +a4* p+\ +2a2* P-\^ = aiix) + asix) , 

h (4-17) 
P-{y) + p_{y)+ai*p-\y + 2a2*p+\y = ai{y) + as{y) • 

Here p±{x) arc the root densities of the two types of strings (4.15), are the corresponding densities 
of 'holes' in these distributions, 

TT 4u^ + 

and * denotes a convolution, i.e. 

/oo 
dx' an{x-x') f{x') . (4.19) 
-oo 

The energy per site of such states in the thermodynamic limit is obtained from (4.8) 

^lirn^ ^-Bjv+,jv_ = -27r ^ j dx p„{x) [ai{u) + a3{u)] . (4.20) 

(7— i 

The eigenstates of lowest energy within this class are obtained by filling the states with negative energy. 
Such states correspond to distributions with Pa{x) = {p'^{x) = 0) for |a;| > A„ {\x\ < A^). They are 
characterized by the total density of type-± strange strings 



L J-A^ 



dxp±{x) . (4.21) 



For the low-lying multiplets identified for small systems in the Appendix B.2.4 we have A^-f + N- w L. 
Due to the symmetry of the equations we can restrict our attention to the case > Within 
the description in terms of string densities p± these states correspond to the choice of = oo: it is 
straightforward to see in this case that ^+^^'- = i independent of A^ . 

We can now proceed and apply standard Bethe ansatz methods to extract the finite size energy 
gaps AE of these low-lying states with respect to the completely filled Fermi seas of type-± strange 
strings (corresponding to A± = oo). Note that to extract the scaling dimensions from this spectrum 
according to the finite size scaling predictions of conformal field theory the gaps have to be measured 
with respect to the true ground state discussed in Section 4.2. This offset between the lowest excitation 
considered in this section and the ground state energy Eq is not accessible to the analytical methods 
applied below and will be determined by numerical solution of the Bethe equations. 

4.3.1 Excitations with A^+ = N_ 

The particular subset of solutions discussed above with = A_ needs to discussed separately. Let 
us first consider the thermodynamic limit. Within the integral equation approach the lowest energy 
state of this type is described by A+ = A- = A and the densities of + type and — type strings are 
found to coincide. Introducing p+{x) = P-{x) = p{x) one is left with a single integral equation instead 
of (4.17). For states with N+ = N_ L/2 such that 

hm ^ = 1 (4.22) 
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this integral equation reads 

/oo 
dx' [2 a2{x - x') + a4{x - x')] p{x') = ai{x) + 03(0;) , (4.23) 
-00 

which is solved by Fourier transformation giving 

p(x) = - — \ . (4.24) 

' 2cosh7ra; ^ ' 

Substituting the solution of the integral equation into (4.20) we then find that the energy for states 
with ^ = I 

lim = _4 . (4.25) 

L — *oo L/ 

All states (of the type we are considering) with = ^ w | are thus degenerate with the singlet 
ground state in the thermodynamic limit. 

Let us now turn to solutions of the Bcthc ansatz equations in the finite volume. Solving these 
equations for finite systems we observe that an important role is played by root configurations in which 
each strange-string of type + becomes "degenerate" with a strange-string configuration of type — . 
After suitable relabelling of the indices this implies yj = Xj for j = 1, . . . , in (4.15). In such a 
configuration the spectral parameters of the two levels of the Bethc ansatz equations (4.6) coincide 

u'f^ = tS^^ , u'^^ = tJ-) . (4.26) 
We will analyze such solutions to the Bethe ansatz equations next. 



4.4 The Takhtajan-Babujian subset 

The results of our numerical solution indicate that such "degenerate" solutions are a subset of an even 
larger class of solutions to the Bethe ansatz equations (4.6) in sectors with same number of roots on 
both levels, i.e. A'' = M. These solutions are obtained by setting all spectral parameters uj on the first 
level equal to the spectral parameters on the second level 7^ , i.e. by requiring that 

Uj = 7j 

The roots then have to satisfy 

/ Uj+i V 
\Uj - i) 

This system of equations is identical with the Bethe equations for the spin-1 antiferromagnetic Takhtajan- 
Babujian (TB) chain with L lattice sites and antiperiodic boundary conditions [56, 57] 

L 

/f^^ = ^(S,-.S,+i-(S,-.S,+i)2) . (4.29) 

Here Sj arc sl{2) spin-1 operators. The degenerate strange string solutions (4.15) with A_|- = and 
(4.26) discussed above become the 2-strings = Uj ± i/2 in the Takhtajan-Babujian-classification. 
We note that the normalization of the Hamitonian (4.29) is such that the gaps in our system will be 
twice the gaps in the TB chain. Alcaraz and Martins have carried out a finite-size scaling analysis 
of the anisotropic TB chain with general twisted boundary conditions [57] and identified the operator 
content of the underlying conformal field theory. We summarize some of their results in Appendix C. 

Boundary conditions must be handled with extreme care. This is especially true in theories with su- 
per group symmetries: since one deals with graded tensor products, periodic and antiperiodic boundary 
conditions for fermions do not necessarily translate into periodic and antiperiodic boundary conditions 
for the spin operators in the chain. 



j = l,...N. 



T-r Uj -Uk+l 

Uj — Uh — i 

k=l J 



(4.27) 
(4.28) 
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4.4.1 Periodic boundary conditions 

The sector we are considering here, = M, is simple in that respect: in that case, since the total spin 
is integer, we are exploring only bosonic states. Therefore, periodic boundary conditions for the 
fcrmions translate into antiperiodic boundary conditions for the spins of the Takhtajan-Babujian chain 
(for a detailed discussion of boundary conditions in chains with supergroup symmetries, see [58]). 

To proceed, we now recall some results about that chain [57] (see also Appendix C). The best 
way to encode the part of the spectrum that comes from the Takhtajan-Babujian chain is to use the 
generating function of levels - the conformal partition function [59] . Let us introduce the quantities 

Anm = ^[n + {m + ^/n)f 

Anm = l[n-{m + $/7r)]' . (4.30) 
o 

We further introduce the partition functions for the Ising model with twisted boundary conditions 



^/sms^Q^O) 



Xo 



2 



+ 



Xl/2 



1 fm . \o. 



2 V \V\ \V\ 



Ising 
Xl/16 



2 _ 1 |^2| 



Z^-"«(0,1) = Zz.,„,(l,0) 2 

z^-"«(i,i) = x5«^"^x{;r+xo^'"'x{;r = UT^ (4-31) 



2 V \v\ \v\ 

The generating function of the levels in the Takhtajan-Babujian chain is then 

C?n = E E ZrT' E g^-""- V24^-A._l/24 (4 32) 

r=o s=o m = r + 2Z 

S + 2Z 

and coincides with the partition function of the SU{2) WZW model at level two [59]. 

We notice that, in the equivalence between the Takhtajan-Babujian spectrum and a subset of our 
sl{2/l) intcgrable spin chain spectrum, the Takhtajan-Babujian model has an even number of sites if 
the sZ(2/l) model has an even number i of 3 (and thus 3) representations. Expression (4.32) is relevant 
to this case only. One can also consider the case of this number being odd. Then the generating function 
of levels reads 

Zif, = E E ^rT' E ,A_-l/24^-A._l/24 (4 33) 

r=0 s=0 m=r+2Z+l,n=s+2Z+l 

In the case $ = tt however, expressions (4.32) and (4.33) can be shown to coincide! 
Let us now extract some conformal weights for the sZ(2/l) chain from these expressions. To do so, 
one has to 

• Double the gaps, because the normalizations of the Hamiltonians are such that the energies of 

the s^(2/l) chain arc twice the energies of the TB chain. 

• Shift the energy differences by a constant ^ (after the rescaling) . This shift in necessary, because 
the finite-size spectrum in the TB chain with antiperiodic boundary conditions is calculated with 
respect to the ground state with L even and periodic boundary conditions [57]. The energy 
difference between the true ground state and the TB ground state for even L and periodic 

2 

boundary conditions is equal to 
For r, s = 0, we get from (4.32) (that is L even) dimensions in the c = theory which read 

, - (m + l)2 1 

h-h = n{m + l); n,m€ 2Z (4.34) 
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Hence we have the possibihty n = 0, m = 0,m = —2, giving rise to h = h = |. Ah other scalar 
operators have h = h = ^ + integer. The lowest non scalar operator corresponds to n = 2, m = — 2 
with /i = 17/8= l/8 + 2,/i= 1/8. 
For r = 1 and s = 

, - (m + If 

h-h = n(m + l); n€ 2Z,m€ 2Z+1 (4.35) 
Here all the dimensions are integers. For r = and s = 1, 

- (m + 1)^ 

h-h = n(m + l); n€ 2Z + 1, me 2Z (4.36) 
with again only integer dimensions. Finally, the case r = s = 1 leads to 
^ - n2 (m+l)2 3 

ft -ft = n(m + 1) + 1; nS 2Z + I, m € 2Z+1 (4.37) 

The case n = 1, to = — 1 corresponds to ft, = 9/8, ft. = 1/8. 

In fact, the generating function of the sl{2/l) dimensions deduced from the TB spectrum can be 
written after a few manipulations as the sum of two simple contributions 

oo 

ZTB subset ^ 2^1/8 "Q (1 + + q2n^ X {q ^ q) (4.38) 

n=l 

and 

oo 

ZTB subset =Y[{1 + g"+l)(l + X {q ^ q) (4.39) 

n=0 

Of course these two contributions boar resemblance to the two Ramond characters discussed in the first 
section. However, it is important to understand that we cannot deduce from these the degeneracies of 
the sZ(2/l) chain. Rather, all we can deduce form the argument is that the the dimensions appearing 
in the expansion of these two objects appear also in the sl{2/l) spectrum, with a multiplicity at least 
equal to the multiplicity it has in this expansion. We denote such a relationship by 

yTB subset ^ lysli^jX) chain 
^1/8 S ^1/8 

^TS subset ^ ^«'(2/l) chain 

It can be checked that this is compatible with the spectrum of the chain containing at least the modulus 
square of the integrable characters, that is 

ZTB subset < |^^^^^^|2 

ZTB subset < |^^JV|2 (4 41) 

We can make some more detailed guess about which states on right are represented by a term on the 
left in (4.40), (4.41). For instance, Z^^'*"''*'^* contains two fields with dimension ft = ft = 1/8, and we 
know that these fields must have B = 0. These fields must also have vanishing SU (2) spin, since they 
are obtained by taking the n number of the Takhtajan-Babujian chain equal to zero. This means that 
their SU{2) spin in the s/(2/l) chain must also vanish, and thus, in the continuum limit, they must be 
the symmetric and antisymmetric combinations of Jq = 1/2, Jq = —1/2 and Jg = —1/2, Jq = 1/2. 
The fields with ft = 1/8 and Jq = Jq = are, in particular, absent in the subspectrum obtained from 
the Takhtajan-Babujian chain. But we know they are there in the full spectrum. 
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Similarly, ^^-Bsubset contains a single field with h = h = 1; once again it must have B = = 0. 
On the right hand side, we have 8x8 fields with such weights, obtained by tensoring the adjoint with 
itself. It is highly likely that the field on the left comes from the part of this tensor product where the 
Jq = 0, Bo = field for right movers is tensored with the Jq = 0, Bq = field for left movers. 

In Tables 1 and 2 we summarize the lowest scaling dimensions (that is, x = h + h, irrespective of 
whether h = h ov not) in the TB subset and their s/(2/l) quantum numbers for L even and L odd. 
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Table 1: Smallest scaling dimensions and s^(2/l) quantum numbers in the TB subset for L odd. 
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Table 2: Smallest scaling dimensions and sl{2/l) quantum numbers in the TB subset for L even. 

We note that these scaling dimensions occur with the same multiplicities as in the odd L case, but the 
sZ(2/l) quantum numbers are different. 

4.4.2 Antiperiodic boundary conditions and the Neveu Schwarz sector 

We can also consider the case of antiperiodic boundary conditions for the spin chain in the sector 
A'' = M. This should correspond to antiperiodic boundary conditions for the fermions, ie the NS 
sector. Note that the sl{2/l) symmetry is then broken. In Appendix B.6 wc summarize some results 
of a numerical solution of the Bethe ansatz equations with antiperiodic boundary conditions in the 
sZ(2/l) chain for small L. 

The same argument we used above in the periodic case shows that the sZ(2/l) chain with antiperiodic 
boundary conditions has a particular sector that is related to the spectrum of the Takhtajan-Babujian 
chain with periodic boundary conditions, i.e. $ = 0. This time, the generating function of the s/(2/l) 
chain levels can be written as the sum of two contributions 

oo 

ZTB subset ^ ^1/4 -Q^-^ ^ ^2n^3 x (g ^ g) , (4.42) 
1 

and 

^TB subset ^ 2q-l/8 ^ (g ^ g) . (4.43) 

Now the correspondence between R and NS weights is as follows 

hNS = /ifl + ^ - Jq'^ , 

hNS = hR + \ + Jo'"" ■ (4-44) 

We thus see that the two fields in Z°^g™ with weight ha = ha ~ 1/8 should get, respectively, 
hNS = Tt-ms = ^1/8, and h^s = Tins = 7/8 = —1/8 + 1. The field in the 1/8 multiplet whose 
dimension becomes Hns = h-NS = 3/8 = —1/8 + 1/2 is the one with Jq'^ = Jq'^ = 0, and it is not 
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present in the Takhtajan-Babujian spectrum. Thus the form of Z'^^j'^ is compatible with the SU(2/1) 
WZW hypothesis, in the NS sector. 

Similarly, the singlet in Z^^ subset j-^g^g h = Jq = and thus acquires, under spectral flow, a 
weight hNS = hNS = 1/4. As for the multiplet with hn = hn = 1, under spectral flow, it acquires 
hNS = hNs = ^/'i- The field with h^ = hji = 1 and Jq'^ = —Jq'^ = 1/2 which acquires under spectral 
flow hNS = ^NS = 3/4 = 1/4+1/2 is absent from the Takhtajan-Babujian part of the spectrum. 

It is also satisfactory to see that, within the Takhtajan-Babujian part of the spectrum, we have 

hNS = ^(n + m)^ = ■^(n + m+ 1)^ - i - ^(n + m) =hR+^ - Jq'^ , 
hNS = ^{n-mf = ^{n-m-lf -j + ^in-m) = hR + j + J^'^ , (4.45) 
and therefore 

Jq^ =^{n + m+l) , 

J^'«=l(n-m-l) . (4.46) 

Observe in particular that Jq'^ -|- Jq'^ = n, the scalar quantum number. All these relations are 
compatible with our WZW conjecture for the continuum limit of the sl{2/l) chain. 

4.5 Strange string excitations with iV+ 7^ N_ 

The integral equations (4.17) for strange string excitations with iV+ 7^ iV_ can be analyzed by Fourier 
transformation. This allows one to express the Fourier transforms of the hole densities p±{oj) in terms 
of the Fourier transforms of the root densities p±{co). Inverting the resulting matrix equation, we 
obtain 



^^^'^^ = 2cosh(a;/2) ' + ^"^"^ ' 



(4.47) 



~ r ^ 1 l + e2|-l ^ ^ pj"! . . 

p-iuj) = , , ^ p_{(jj) H 5— p+Ujj) . 

^ ' 2cosh(w/2) 4sinh2c<; ^ ' 2ii\uh^ u ^ 

Here we have defined the Fourier transform by f{uj) = J dxe^'^^f{x) and a„(w) = e""!"!/^. The 
right-hand side of (4.47) is most conveniently written in matrix form using 

The matrix R is the Fourier transform of the resolvent R= {1 — K)~^ of the (matrix) integral operator 
in our original equations (4.17) with kernel 

^ \ 2a2{x-y) a4{x-y) J ^ ' 

Starting from (4.17) (or, cquivalcntly, (4.47)) and (4.20) wc can analyze the finite size spectrum of these 
excitations by standard methods. Cases similarly to the one considered here with complete symmetry 
between the excitations in the + type and — type sectors have been studied in Refs. [60, 61]. There 
the following general formula for the finite size spectral gaps AE of low lying excitations over a filled 
Fermi sea with several components in a relativistic invariant model (i.e. with a single Fermi velocity) 
was derived. In our case it takes the form 

—AE = -(ANVR-^(0)AN + D^R{0)D . (4.50) 
2nv 4^ ' ^ ' 



Here u = tt is the characteristic velocity of the system and is calculated in Appendix B.5, AA'' is a 
MA^V 
\AN_ 



2-component vector (^^j'^ ] where ANc^ is the change in the number of roots of type a compared to 
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the ground state. Similarly, D is the 2-component vector ^^^^ where is the number of solutions 

of type (T "backscattered" from the left to the right of the Fermi sea. 

In our case there is a slight additional subtlety here coming from the fact that the ground state is 
made of complex solutions (4.15): in such a situation it is known that - even in the usual cases such 
as sl{2) - on top of the Gaussian spectrum, deduced from our general formula (4.50), the dimensions 
of other operators may appear, as well as selection rules. At the same time it has been observed in 
these cases that most of the basic qualitative conclusions drawn from an analysis based on the string 
picture are valid and provide the basis for a quantitative description when amended for example by 
results from a numerical computation of these gaps (see e.g. [57, 62, 63]). 

Applying (4.50) blindly to our case requires a bit of care as the matrix R has divergent elements. 
This rather unusual fc^aturc can be traced back to the divergence of the Fourier transform of the kernel 
describing the scattering between complexes of type + and — for — iV_. We can regularize these 
elements by giving oj a small value 

Prom this it follows that 

—AE = -{AN+ + AN_f + X {AN+ - AN_f + -{D+ + D_f + oox{D+- D_f (4.52) 

2'KV 2 8 

A divergence of the kernels similar to the one leading to (4.52) here has been encoutered in the study of 
the osp{2/2) chain (that is, based on the four dimensional typical representation [0, 1/2]) and the study 
of the antiferromagnetic Potts model. In both cases, the divergence of the kernels (in Fourier space) 
at w = has been interpreted [31] as a manifestation of a continuous spectrum - which is equivalent 
to the infinite degeneracy of the finite size gaps between states with N+ + N- = const, in (4.52). 

We emphasize that this result is obtained for 0(1) deviations in N± compared to the state with 
equal numbers of the + and — type strings (4.23). In other words, the result holds for finite values of 
AA^, D. 

In fact, the state with the largest possible asymmetry between these distributions, i.e. Nj^ = L, 
N- = 0, can be analyzed in the thermodynamic limit giving 

, , 1 cosh^a; , , , , 

for the densities. The energy of this state is Sjv+,o = — 2L(l+ln2) and the state is clearly highly excited. 
In summary, we expect on the basis of the above analysis that there is a macroscopic degeneracy of 
the finite size spectrum identified in the Takhtajan-Babujian subset. 

As we have previously mentioned, these results ought to be taken cum grano salis as they have been 
obtained within the framework of the string hypothesis. The latter neglects the finite size deviations 
e± in (4.15), which also contribute to order 0{L~^) to the gaps. Furthermore, in order to extract 
conformal dimensions from the finite size spectrum we need to determine the difference between the 
groimd state energy Eq (see Sect. 4.2) to that of the state AA^± = = D± to order 1/L which is 
beyond the scope of the analytical methods applied here. In light of these complications and in order 
to verify the observation from (4.52) that each gap is infinitely degenerated (or equivalently that the 
spectrum has a continuous component added to the discrete part identified via the Takhtajan-Babujian 
mapping), we have carried out an extensive numerical analysis of the Bethe ansatz equations (4.6) for 
large, finite L. Due to the presence of strong logarithmic corrections (i.e. corrections of order 1/(L In L) 
to (4.50) we had to study systems with L up to 5.000, i.e. 10.000 lattice sites. 

We summarize the results of this finite-size scaling analysis for several low-lying [0, 1] octet states 
+ N- = L — 1) as well as the octet and the indecomposable [0,-1/2,1/2,0] representations 
(A'^-i- = = L/2) in the Takhtajan-Babujian sector in Appendix B.4. The scaling of their energies 
expected from conformal invariance is 

^[E{N+,N_)-Eo] >h + h = x, (4.54) 

ZTTV I/— >(X) 
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where x scaling dimension. We plot the left-hand side of equation (4.54) as a function of the lattice 
length for several different excited states in Fig. 5. 




, jl , \ , \ , L 

0.2 0.4 0.6 

mog{L) 



Figure 5: Corrections to scaling of the energies of the low- lying states studied in Section B.4. The 
dashed lines indicate the rational function extrapolation L ^ co. 

Within rational function extrapolation of our numerical data wc find that the energies of the statc;s 
considered all become degenerate to order 1/L in the thermodynamic limit. The scaling dimension 
obtained from this extrapolation of (4.54) is 




(4.55) 



As has been discussed above, we have to extend the analysis of the finite size gaps to the case of 

antipcriodic boundary conditions for a proper identification fo the continuum limit. For this we have 
studied states build from the strings (4.15) in some detail which again are the relevant configurations 
for the lowest excitations. The details of the numerical analysis are given in the appendix, here we 
just present the main results which give some evidence that again there are are many states degenerate 
with the ground state (in the Takhtajan-Babujian sector) in the thermodynamic (L oo) limit. All 
of these states have A^_|_ -I- A''_ = L, but different values of AiV = A'^^ — N_. Their scaling dimensions 
extrapolate to x = — 1/4 relative to the ground state of the chain with periodic boundary conditions. 
The corrections to this scaling for some states are shown in Fig. 6. 

4.6 Summary 

As a summary, we represent on figure 7 the lowest gaps we have observed in the R and NS sectors. 
Values in parenthes is are expected from the study of integrable characters of the WZW model, but 
we have not been able to investigate them numerically. 

5 Higher level 

It is interesting to consider the same problem with larger representations - more precisely, the fully 
supersymmetrized atypical representations of type [±j, j] [ s = 2j boxes in the Young diagram repre- 
sentation). In this case, the Bethe equations allow for a sector reproducing twice the gaps of the spin 
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Figure 6: Finite-size scaling of low-lying energy levels for the system with antiperiodic boundary 
conditions. Eq is the ground state in the sector with periodic boundary conditions (see Section 4.2). 



s Takhtajan-Babujian chain. For the latter, the gaps are obtained once again by combining the gaps 
from a free boson and disorder operators from a Z2s parafermionic theory. The gaps from the free 
boson sector read now 



1 2 
A„m = — [n + {m + s^/n)] 

1 2 

A„„ = — [n - (m -t- s$/7r)] 



(5.1) 



The sl{2/l) symmetric case (R sector) corresponds to $ = tt. The leading gaps with respect to the 
c = theory read then 



h + h = 



(m + s)^ s r(2s — r) 

2s 2(7+1) 2s{s + l) 



h = h 

with the constraint that m = r mod 2s. The choice m = r — 2s gives then 

- (s - r)2 



(5.2) 



(5.3) 



This obviously reproduces the conformal weights in the SU{2/1) WZW model for representations [0, j] 
with j = 0, 1/2, . . . , ^ at level s. 

Meanwhile in the Neveu Schwarz sector, Cg// = ^fj, and for the weights we get the other formula 



h = h 



m 



4(s + l) 



+ 



(m — s + r)(m + s — r) 
4s(s + l) 



The case m = r — 2s now gives 



h = h 



(2s-r-H)^ _ 1 
4(s-t- 1) 4 



(5.4) 



(5.5) 



28 



teS/B 
h=1 



h=ly8 



(h=3/8)- . 
h=1/4 - ■ 



c=0-^ h=-1^ 

c, = 3± 



NS 



Figure 7: The spectra from the lattice model in the R and NS sector, in agreement with the WZW 
prediction. The thick lines represent a continuum of critical exponents (extending to infinity). Values 
in parenthesis were expected but beyond our reach. 



This can be written as 

matching the formula for the spectral flow in the SU{2/1) level s WZW theory. Hence central charges 
and TB subset match the WZW prediction in this case as well, and it is likely that a more thorough 
study would find similar agreement for the continuum parts of the spectrum too. 



6 Conclusions and speculations 

To conclude, from the information at our disposal using the Takhtajan-Babujian part of the spectrum, 
analytical study of the Bethe ansatz, and extensive numerical calculations, we have obtained consistency 
with the hypothesis that the continuum limit is the SU{2/\) level one WZW model. This consistency 
includes the set of all dimensions in the R sector, and only a subset of dimensions in the NS sector. It 
is highly indicative that our conjecture is correct that the effective central charge in the NS sector is 
c = 3. As for the part of the NS spectrum whose weights are of the form Jins = —\ + \ + integers 
hNS = i + i + integers, we have argued they are situated outside of the Takhtajan-Babujian part 
of the spectrum, and thus we have nothing to say about them based on this chain of arguments. It 
would of course be a crucial further test to extract these exponents from the Bethe ansatz solution 
(this might be difficult as there are sitting inside the continuum). We note that it would be extremely 
useful to compare more thoroughly than we have done the group content of the lattice and field theory 
models: this however may be a very non trivial exercise, in particular because fine structures of the 
representations (such as indecomposability) may not evolve straightforwardly in the continuum limit 
process [65]. 

One of the main difficulties we have encountered is the the lack of understanding of the WZW 

theory itself. In particular, the free field representation as well as algebraic arguments point to the 
existence of a spectrum of dimensions unbounded from below, and maybe continuous. While this is 
somewhat expected in view of the unboundedness of the action in the field theory, it is not totally 
clear how this should be related to the spectrum and partition function of a lattice model. We have 
suggested that some sort of analytical continuation of the characters is at work, folding back the 
spectrum of arbitrarily large negative dimensions into a continuum of positive dimensions. This is 
strongly supported by calculations involving modular transformations of characters, and gives rise to 
predictions in remarkable agreement with lattice calculations. A point that remains unclear is whether 
this continuation is a formal operation, or has the meaning of changing the target space to a Riemannian 
symmetric superspace as argued in [33]. Clearly more work is needed to clarify these questions, which 
we believe are crucial to our understanding of CFT description of the transition between plateaux. 
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Finally, we come back to our original question concerning the universality of3(g)3ig)3ig)3... sl{2/l) 
spin chains. It is easy to check that the analysis in [16] for the "Heisenberg" model (with S.S nearest 
neighbour coupling only) in the R sector, gives a totally different spectrum. First, multiplicities are all 
finite. Then, the following scalar operators 

, - (3P+ 1)2-1 , - 4M2-1 

h = h=^ , h = h= — — (6.1) 

24 24 ^ ^ 

appear, as well as many other non scalar operators whose chiral weights generically read (modulo 
integers) 

Except for the lowest weights, there is no indication that any of these weights appear in the spectrum 
of the integrable s^(2/l) chain. Moreover, imposing antiperiodic boundary conditions for the fermions 
(ie the NS sector) gives results which are totally different from the ones of the integrable spin chain; 
in particular, instead of being three, the effective central charge is an irrational number c = 1 + 
_^ In^ ^'^^ PS 1.84! In fact, pretty much the only features common to the two systems are the 
vanishing of the central charge in the R sector, and the lowest gap h = h = ^ - extremely ubiquitous 
features of superalgebra field theories. 

Since the integrable model requires fine tuning of nearest neighbour and next to nearest neighbour 
couplings, we expect that it describes the least stable of the two fixed points, as evidenced by the 
higher value of its effective central charge in the NS sector. 

It seems desirable to investigate the complete phase diagram of the model, a task that could be 
made easier by using the loop representation. We hope to get back to this question later. A property 
that could prove useful in this analysis is that the ground state of the Heisenberg chain, like the ground 
state of the integrable spin chain, is a true singlet, with trivial value of the energy. It does not seem 
to be the same state however. 
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A Some algebraic considerations 

Numerical study indicates that there is a unique singlet in the product (3 (g) 3)®-^: for instance, the 
tensor product (3 55 3)*^^ decomposes as 

(3 ^ 3)2 = [0, 0] + 4[0, 1] + [0, 2] + [1/2, 3/2] + [-1/2, 3/2] + [0, -1/2, 1/2, 0] (A.l) 

One has to be careful in defining what one means by singlet in this problem, as there are invariant states 
(that is, states annihilated by all the generators of the algebra) within some of the indecomposable 
representations too. We will call 'true singlet' a representation isomorphic to the trivial representation, 
and which is not a quotient of a larger indecomposable representation. When decomposing (3 ® 3)®^, 
it appears as an invariant state which does not lie in the image of the raising or lowering generators of 
the algebra. 

It is known from representation theory of sl{2/l) [24] that in the tensor product we are interested 
in, there can only appear typical representations, simple (irreducible) atypical of the type [±J, j] and 
non simple (indecomposable) atypical, representations. The latter have vanishing superdimension, as 
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do the typical representations. It is obvious that the superdimension of our tensor product is one but 
one could imagine that this one does not come from the existence of a single [0, 0] representation. For 
instance, it could arise as the result of a cancellation between several [0, 0] reps (each of superdimension 
one, as it is easy to see that all [0, 0] states have to be bosonic in our system) and several simple atypical 
representations with negative superdimension. It is easy to exclude this possibility if one recognizes 
that the typical representations as well as the indecomposable atypical appearing in our tensor product 
are projective, while simple (irreducible) atypical are not (recall here [28] that a representation is 
projective if it cannot appear as a quotient of a bigger, indecomposable block). In particular, [0,0] is 
not projective. A well known result of algebra [66, 28, 67] says that 



But consider 3 3 = [0, 0] + [0, 1]. The adjoint appearing here is typical and thus projective. Hence, 
tensor products [0, 1]®-^ decompose onto projoctives only, that is typical representations and non simple 
atypical. No simple atypical ever can epppear in these tensor products. Thus, [0, 0] can appear only 
once, as the result of the tensor product of [0, 0] with itself, which is our initial claim. 

Note that the result could also be proven by using the fact that the tensor product of any in- 
decomposable representation (either irreducible like a simple typical, or reducible like a non simple 
atypical) of vanishing superdimension decomposes as a sum of representations each with vanishing 
superdimension themselves [68]. 

Notice that a similar result is known to take place in the SU{2)q product of an even number of 
fundamental (two dimensional) representations when = —1, so the spin j = 1 representation has 
vanishing q-dimension. In this case, there is a unique singlet representation, which coincides with the 
space Ker S+ /Im S+ [69, 70]. 

Of course, if we get back for instance to (3 (g) 3)®^, there are two invariant (that is, annihilated by 
all the generators in the algebra) vectors, one in [0, 0] and one in the indecomposable. The vector in 
[0,0] is a 'true singlet', which means simply that it is not in the image of the generators. Obviously, a 
generic combination of the two invariant vectors will still have this property. This means that the way 
[0, 0] will be expressed in terms of the basis vectors in 3 and 3 will in fact depend on the Hamiltonian 
(one could imagine selecting the 'good true singlet' by requiring that it be orthonormal to the image of 
the generators in the subspace of vanishing U{1) charge. However this will not work, as the invariant 
state in the indecomposable has zero norm square) . Remarkable expressions for such states have been 
obtained recently (albeit in a different language) in the context of the 'Razumov Stroganov conjecture' 
[71] . We leave this point to further study. 

B Details of the analysis of the lattice model 

B.l String hypothesis 

The standard way to classify the solutions of (4.6) is to consider configurations of spectral parameters 
that sit on poles of the bare scattering phase shifts (the right-hand-sides in (4.6)). A simple calculation 
yields the following types of "strings" 

(1) "Reals": 

unpaired, purely real spectral parameters uj and 7/3. 

(2) "Wide strings": 

"type-1": composits containing n — 1 7's and n u's (n > 1) 



projective (g) anything = projective 



(A.2) 



7, 



u 



S^f = ^,^1) + i + 1 _ 2fc) , k = l...n 
'af = u^^''^+i{n-2j), i = l...n-l, ni"'i)eIR 



(B.l) 
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"type-2": composits containing n 7's and n — 1 u's (n > 1) 



(n,2) 
(n,2) 



7^"'^^ + i (n + 1 - 2fc) , = 1 , 
7^2) + J („ „ 2j) , j = l...n 



1, 7i"''^G]R 



(B.2) 



Figure 8: Wide strings of lengths three and two respectively. The circles/triangles denote the positions 
of the u's/7's involved in the string. 



(3) "Strange strings" : 

composits containing n 7's and n it's (n > 1) 



or 



(" -) 
(" -) 



7a,fc 



ul^'> +i{n+l 
-i{n+l 



u 



i(n + 1 



2fc 
2fc 

2k 
2k 



k = l,... 

fc= 1,... 



(B.3) 



(B.4) 



We note that solutions of this type of solution is quite different from usual string solutions in that 
the set of roots on one level of the Bethe equations is not invariant under complex conjugation. 
The other solutons discussed above (reals and wide strings) are invariant under this operation. 
This is somewhat similar to what was found recently in [55] . It is important to note that although 
"strange strings" are not invariant under complex conjugation, the corresponding energy (4.8) in 
the case A = is still real. 

(4) "Narrow strings" : 

composits containing n 7's and n w's (n > 1) 



ur-^' + -{n+l-2k), k = l. 



{n,n) 



4"'"^ + 2("+l-2j) 



,7 = 1 



. n 



(B.5) 



Narrow strings may be thought of as special cases of strange strings or wide strings in the following 
sense. 
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(a) 



(b) 



Figure 9: (a) The two types oi n — 1 strange strings; (b) A "+" type strange string with n = 2. 



Combining a "+" -strange string of length n and centre w*^"-' with a a "-"-strange string of 
length n and centre m^"-* we obtain a narrow string of even length 2n. This is shown for 
n = 1 in Fig. 10(a). 

Combining a type-1 wide string of length n and centre u^"'^-' with a type-2 wide string of 
length n and centre w*^"'^^ we obtain a narrow string of length 2n — 1. This is shown for the 
case n = 2 in Fig. 10(b). 



ik Im 



(a) 



Re 



(b) 



Im 



Re 



Figure 10: (a) Combining a pair of and "-" strange strings of length 1 gives a narrow string of 
length 2; (b) Combining a type-1 wide string of length 2 with a type-2 wide string of length 2 gives a 
narrow string of length 3; 



It is clear from our discussion that narrow strings are not "fundamental" but are merely degen- 
erate cases of strange string solutions of the Bethe ansatz equations. 

We have verified explicitly by numerical solution of the Bethe ansatz equation (4.6), mainly for the 
case A = 0, that solutions of the above types exist at least for small string-lengths n. On the other 
hand it is entirely possible that the above classification is not complete. 

However, it is important to note that if 

(a) our classification is complete and 

(b) in the thermodynamic limit solutions of the Bethe ansatz equations (4.6) with the same values 
of per site and B per site as the ground state tend towards ideal string solutions (as is often 
the case in Bethe ansatz solvable models), 
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then all energies are real for A = 0. Given the string hypothesis (B.1)-(B.4) the usual way of proceeding 
in a Bethe ansatz solvable model is to construct the thermodynamic Bethe ansatz (TBA) equations at 
finite temperature and then use them to determine the distribution of roots in the zero temperature 
ground state. This applicability of this procedure to our case is questionable as it is not clear whether 
the minimization of the free energy functional is a meaningful procedure for the case of a non-hermitian 
Hamiltonian. Disregarding this issue one finds that the standard "integer method" of classifying 
solutions of the Bethe ansatz equations does not seem to work. In particular, we find that in some 
cases the numbers of solutions of a particular type for fixed L, N and M are found to differ from 
what the integer method predicts. Wc don't have a good explanation for this fact at present and leave 
this issue for future studies. Irrespective of this, the Bethe ansatz does not provide a complete set of 
states in any case. This follows from the peculiar representation theory of superalgebras. It was shown 
in [38] that for the model at hand all cigcnstatcs obtained from the Bethe ansatz are highest weight 
states of sl{2/l). For integrable models exhibiting a classical Lie- algebra symmetry a complete set 
of states can be obtained from these highest-weight states by acting with all possible combinations of 
lowering operators [72] . In the present case there arc states that belong to reducible but indecomposable 
respresentations, for which the highest weight state is not a cyclic one [24]. As a consequence, acting 
on the Bethe ansatz states with lowering operators is not sufficient to obtain a complete set of states. 
This complication is absent in other integrable, supcrsymmetric models, for which complete sets of 
states have been obtained by the aforementioned procedure [49, 73, 74, 75, 76]. 

Given these complications we will follow a different, less ambitious approach: by studying finite 
systems we will try to identify the ground state and some low- lying excited states of the Hamiltonian 
(4.5) in terms of distributions of roots of the Bethe ansatz equations (4.6). 

B.2 Small Finite Systems 

For L = 1, 2 the transfer matrix T:i{v) can be diagonalized directly. 



B.2.1 L=l 

The tensor product 3 (8) 3 is fully reducible and decomposes into a singlet [0] and the typical octet 

[&,i] = [o,i]. 

3 ® 3 = [0, 0] [0, 1] (B.6) 
The corresponding eigenvalues of tz{v) are readily obtained 

2 2 



A8(^;) 



Ai(v) = l-— — -- + 



A-12 A-1 2 (^-^^ 



\+l V A+lu-l-A 
Similarly we find for the eigenvalues of 73(11) 

2 2 



V V-l + \ „x 

, ^ ^ A + 12,A+1 2 (^-^^ 

^^(^^ = ^-a^^ + a^^;3tta- 

Obviously, A8(f), A8(u) are the eigenvalues on the pscudo vacuum (M = = N) while hi{v), Ai(w) 
can be obtained from (4.7) in the sector with {M,N) = (1, 1) for m = A, 71 = which clearly solve 
the BAE (4.6). We note that the singlet eigenvalue satisfies the functional equation 

A.MA.(,-l-A).( '°^;;_Y:-;-^' )\ (B.) 

Furthermore, for A = 0, we find Ai(w)|a=o = {v — 2){v + l)/v{v — 1) = Ai(w)|;^=o- The eigenvalues of 
the Hamiltonian are 

^1 = , ^« = A^- (B-IO) 
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B.2.2 L = 2 



The tensor product of representations appearing in the L = 2 [A site) system can be decomposed as 

(3 ® If = [0, 0] + 4[0, 1] + [0, 2] + [1/2, 3/2] + [-1/2, 3/2] + [0, -1/2, 1/2, 0] (B.ll) 

Here, the atypical representation [0, —1/2, 1/2, 0] is an eight dimensional reducible but indecomposable 
representation which is the semi direct sum of [1/2,1/2], [—1/2,1/2], and two [0,0] representations 
[24]) (see section 2 of the main text). Note that this block is the same indecomposable block as the 
one which appears in the product [0, 1/2] (g) [0, 1/2]. 

Diagonalization of the transfer matrix yields the eigenvalues of the 4-site system. Below we present 
their values together with the labels of the sectors (M, N) in which they appear in the BA and the 
corresponding irreducible representation of the algebra sZ(2/l) in the decomposition (B.ll): 



(M, N) 


[B,S] 


eigenvalue of t3{u),t-^{u) 


Bethe roots 


(0,0) 


[0,2] 


AM = 1 * + 4 + f X + 1-^^ 

A /'„,^ _ 1 4 1 4 1 4 1 4 
AielWj - 1 „ + „2 + „_i+A 1 („-l+A)2 


VcLCUum 


(0,1) 




Ai+^(t;)-l 4 4+ fx + , / 

-1-2 \ ^ i; V'^ 17— 1 — A (U— 1 — A)^ 

Ai2 (w) - 1 .„ + ^ + „_i+A («-l+A)a 


Ml = 


(1,0) 


[_i 31 
L 2' 2j 


k\:\v) = 1 4 + 4 + 4 4 
A-12 (w) - 1 „ „2 + „_i+A ' (?;-l+A)2 


71 = -^'^ 


(1,1) 


[0,1]' 




Mi = -i[A±VA2-l] 
7i = ui + iA 


(1,1) 


[0,1]' 




Zi,= i[^±^3^^+3] 

7i = — iA — Ml 


(1,2) 


[0,-iio] 


- 1 ^ ^ + (v-\-X)'^ 
KJV) = \ 4 4+ f.ux ^-iOT 


ui,2 = -i[A± Vl + A^] 
71 = -i\ 


(2,1) 


[0,-iioj 


A,(«) = 1 ^ 1, + * r i%^2 

A2:(w) — 1 „ + y^i^x (t;-l + A)2 


Ml = 

71,2 = T«Vl + A2 


(2,2) 


[0,-iioj 


V / f ' v — \ — X {v — X — aY 

K[v) = \ 4 4+ f^x ^^XW 

^ f i;^ f — 1+A (u— 


■"1,2 = J_, J_ 

71 75 


(2,2) 


[0] 


Ai(^)=(^(i^ir-yT 
Ai(^)=(^^;gir;Ar)' 





Table 3: Spectrum of the transfer matrices for i = 2. 



Here the eigenvalues corresponding to the four (typical) octets [0, 1] are 



(1,2),, . /A-1 4 A 4 /A-1 



1 + At; VA + 1v2 l + At;-l-A VA + l(t;-l-A) 



^2 



.(1,2),, , A 4 /A+1 4 A 4 ^ /A + 1 4 

^« = ^-a^^^Va^^ + a3T^;3tta^Va^ (^;-i + a)^ 

,(3,4),x X 1 A2 + A + 4±2V)?T3 4 2±\/)?T3 4 
^« = ' OTAp 1 + A ^ 

A^ + A + 4±2\/A2T3 4 2 ± ^/FT3 4 

(1 + A)2 v-l-A 1 + A (w-l-A)2 ^ ■ ^ 

Af^)(A,.) = Af^)(-A,.). 
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Let us consider the [0, —5, 5, 0] multiplet in more detail. There are altogether three highest-weight 
states, i.e. states annihilated by the sl{2/l) raising generators Their quantum numbers are B = ±5, 

= \ and B = ,P = Q respectively. From the tabic above wc sec that the Bcthc ansatz yields three 
states in this multiplet, which by [38] are highest weight states. The Bethe ansatz states with (M, A'') = 
(1, 2) and (M, N) = (2, 1) are readily identified as, B = = \, B = -J^ = -\ highest-weight states. 
The solution to the Bcthc ansatz equations with (M, N) — (2, 2), leading to the eigenvalue A^^, is very 
peculiar to say the least: it contains a free parameter! For any value of 72 the roots given above satisfy 
the BAE (4.6) and give the eigenvalue Kx{v) of the transfer matrix! We identify this last highest- weight 
state with the invariant singlet state of [0, — ^, ^, 0]. The second singlet state in the multiplet is a cyclic 
state, i.e. all states in the representation can be obtained from it by acting with symmetry generators. 
However, it appears that it is not possible to obtain this state from the Bethe ansatz! In other words, 
the Bcthc ansatz is manifestly incomplete due to the peculiar representation theory of superalgcbras. 

For the eigenvalue Ai(i') corresponding to the unique regular singlet in the tensor product there 
exists no non-degenerate solution of the BAE (4.6). Applying a twist in the boundary conditions and 
studying the evolution of BA roots as this twist goes to zero we were able to verify that the singlet is 
in fact described by the singular solution = A, 7^ = which clearly gives Ki{v) from the general 
form (4.7). This finding coincides with the observation in [54] about the distribution of BA-roots for 
the ground state of the quantum transfer matrix for the t J model in the thermodynamic limit. For 
finite L the authors of [54] find this degeneracy to be lifted which is a consequence of the fact that 
their system is subject to a phase shift when compared to ours (see Appendix D). Note that this phase 
shift also breaks the supersymmctry of the system. Furthermore it is not obvious that the two states 
actually are mapped onto one another when the phase shift is changed adiabatically: the distribution 
of roots of the function T) in [54] seems to be very difi'erent. 

We note further that Ai(w) is just the square of the eigenvalue obtained for the singlet in the 2-site 
system (B.7). Hence, Ai(?;) is again the solution of the functional equation (B.9). This observation 
lead us to formulate the following conjecture for the ground state eigenvalue of the transfer matrix 



v{v-l- A) 

For the L = 2 system the eigenvalues of the Hamiltonian are (4.8): 



r(3,4) _ -2A8 - BA^ + hQX'^ + 40A2 + 42 ± V3TA2 (-SA^ + 28A'' + 2Q\^ + 24) 



(A2± 1)3(2 ±V3TA2)2 



For A = we find that E, 



(3,4) 



= —6, i.e. the eigenvalues are equally spaced integers. We note that 

the conjecture (B.13) implies that the ground state Eq = 4,L/(}? — 1) is strictly extensive without any 
finite size corrections. For a conformally invariant system this implies that the central charge of the 
underlying Virasoro algebra is c = — as expected for the class of models considered here. 

B.2.3 L=3,4,5 

By numerically diagonalizing the transfer matrix for system sizes up to L = 5 we find that in general 
there is exactly one eigenstate with an eigenvalue given by 

This is in agreement with the conjecture (B.13). For A = this state is the ground state of the 
Hamiltonian with energy 

Eo = -4X . (B.15) 

From now on we set A = 0. The first few low-lying multiplets for L = 3 are shown in Table 4. We 
see that the first excited state is an octet with energy —11.1231. Then there are two degenerate octets 
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[B,S] 


E 


BA roots 


P/27r 


[0,1] 


-11.1231 


ui = 71 = 0.530i, W2 = 72 = -0.530i 





[0,1] 


-9.60555 


^1=7^ = -0.254 - 0.608i, U2 = 73 = 0.254 - 0.608i 





[0,1] 

L ' J 


-9.60555 


ui = 7i* = -0.254 -h 0.608i, ^2 = 7? = 0.254 + 0.608i 





[0, T7,-T^,01 


-9.00000 


71 = -7* = -0.144 - 0.382Z, ui = -1.427i 
U2 = -U3 = -0.163 + 0.713i 





[0,-i i 01 


-9.00000 


71 = -7* = -0.144 + 0.382i , ui = 1.427i 
U2 = -ul = -0.163 - 0.713i 





[0,1] 


-6.5 


ui = 7i* = -0.238 - 0.547i, ^2 = 7! = 1.104 + 0.175?: 


-1/3 


[0,1] 


-6.5 


ui = 7* = 0.238 + 0.547Z, = 7^ = -1.104 - 0.175i 


1/3 


[0,1] 


-6.5 


ui=jt = -0.238 -t- 0.547i, U2 = 7^ = 1.104 - 0.175i 


-1/3 


[0,1] 


-6.5 


Ml = 7* = 0.238 - 0.547i, ^2 = 7! = -1.104 -|- 0.175i 


1/3 


[0,1] 


-6.23607 


ui = 71 = -0.535 - 0.309f, U2 = J2 = -0.535 -1- 0.309i 


-1/3 


[0,1] 


-6.23607 


ui = 71 = 0.535 - 0.309i, U2 = 72 = 0.535 + 0.309i 


1/3 



Table 4: Low-lying multiplets for L = 3. 



[B,S] 


E 


BA roots 


[0,1] 


-15.0484 


ui = 7* 0.227 + 0.554i, ^2=7^ = -0.227 + 0.554i 
Us = 7| = -0.50948i 


[0,1] 


-15.0484 


ui = 7* = 0.227 - 0.554i, = 7l = -0.227 - 0.554i 
U3 = 7| = 0.509i 


[0,-^,^,0] 


-14.1677 


ui = 71 = -0.051 -1- 0.521i, U2 = 72 = 0.871 + 0.700i 
U3 = 73 = -0.051 - 0.521i, U4 = 74 = 0.871 - 0.700i 


[0,1] 


-13.0749 


Ml = 7^ = -0.426 - 0.654Z, U2 = 72 = 0.426 - 0.654i 
U3 = 7| = -0.596i 


[0,1] 


-13.0749 


ui = 7i* = -0.426 -1- 0.654i, ^2 = 7* = 0.426 + 0.654i 
M3 = 7* = 0.596i 



Table 5: Low-lying multiplets for L = A. 



with energy E = —9.60555. Finally there are two degenerate [0,—^, ^,0] multiplets. We only list one 
Bethe ansatz solution for each of these representations. Like in the two-site case L = 2, the invariant 
singlets of [0, — ^, ^, 0] multiplets correspond to one-parameter families of solutions to the Bethe ansatz 
equations. 

In Table 5 we list the low-lying sl{2/l) multiplets for L = 4. The first excited states are two octets 

with energy E = —15.0484. Next there is a sing Ic [0, -i,i,0] multiplct with energy E =^ —14.1677. 
Clearly the degeneracies for L = 4 are different from L = 3. This suggests that there might be a 
difference in structure of the low- lying excitation spectra for odd and even L. 

Table 6 summarizes the lowlying states for L = 5. The first excited multiplct is an octet, then 
there are two degenerate octets, followed by two degenerate [0,—^, ^,0] multiplets. The structure of 
degeneracies coincides with the L = 3 case. 

B.2.4 Root configurations for a class of low-lying states 

Our results for small systems indicate that the lowest-lying excitations of the model are in the octet 
sector, i.e. belong to [B, S] = [0, 1], or form indecomposable [0, — ^, ^, 0] multiplets. Within the Bethe 
ansatz for a system of length L the octets are parametrized by i — 1 roots Uj and 7q, while the 
indecomposables are found for example in the sector B = = with L roots on both levels. Within 
our classification of solutions by means of the string hypothesis from Section B.l the corresponding 
roots form certain sequences of strange strings of length n = 1 (or their degenerate realization as "2 — 2" 
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[B,S] 


E 


BA roots 


[0,11 


-19.4230 


ui = 71 = -1-0.200 -1- 0.523i, U2 = J2 = -1-0.200 - 0.523i 
M3 = 73 = -0.200 + 0.523z, U4 = 74 = -0.200 - 0.523i 


[0,1] 


-18.66651 


ui = 7* = -h0.390 - 0.593i, ^3=7! = -0.559i 
U2 = 7I = -0.390 - 0.593i, U4 = 7* = +0.503i 


[0,11 


-18.66651 


ui = 7i = -0.390 + 0.593i, U3 = 73 = -h0.559? 
tt2 = 7I = +0.390 + 0.593i, U4 = 7| = -0.503i 


[0,-^,^0] 


-18.28188 


(A) 

ui = 7i* = -1.105 + 0.671i, W4 = 74 = -0.835 - 0.743i 
M2 = 72* = -0.113 - 0.543«, U5 = 7| = 0.02364 + 0.50539i 
M3 = 7| = 0.263 - 0.551i 


[0,-^,^,0] 


-18.28188 


(A) with Uj — > —"^ii la — ^ ~la 



Table 6: Low- lying multiplets for i = 5. 



narrow strings of length (B.5)). We denote the number of these "+" and "— " type strange strings by 
and iV_ respectively. We constrain ourselves to the case AL|_ > A''_; states with A''_ > are 
obtained by exchanging m's and 7's and are by construction degenerate. We denote the energy levels 

of the states by Emj^^N-- 

4.4.1 Even L: 

• A low-energy state with B = Q — j'' belonging to a [0, — i, i, 0] multiplet is obtained by choosing 
N+ = N- = -J and considering a realization of this configuration as a collection of "2 — 2" narrow 
strings. This implies that after a suitable relabelling of roots we have Uj = jj, for j = 1,. . . ,L. 
Using this degeneracy the Bethe ansatz equations (4.6) can be rewritten as 

= n^^^^^, j = l,...,N, (B.16) 

with N = L. These equations are identical to the Bethe ansatz equations of the SU(2) symmetric 
Takhtajan-Babujian spin-1 chain [56] for a lattice of L sites and antiperiodic boundary conditions. 
The operator content of an anisotropic v(^rsion of this model (with twisted rather than antiperiodic 
boundary conditions) was determined by Alcaraz and Martins [57]. We show in Appendix C how 
our results relate to theirs. 

The lowest lying excitation has = ^ and A''_ = ^ — 1 and is the highest weight state of [0, 1] 
multiplet. The distribution of strings is symmetric around the imaginary axis, i.e. 

{uj\j = 1, . . . L - 1} = {-u;\j = 1, ... i - 1}, (B.17) 

and similarly for {"fjlj = 1, ... i — 1}. The root distribution for this state for L = 10 is shown in 
Fig. 11 (a). 

Other [0, 1] octet states are obtained by taking Af+ = ^ + n and = ^ — 1 — n, n = 1,2,..., 



i.e. by replacing n "— " type strange strings by "+" type strange strings. The corresponding 
root distributions are shown in Fig.ll(b)-(f). We note that the lowest energy state for given N± 
corresponds to a root distribution that is symmetric aroimd the imaginary axis. Higher energy 
states correspond to asymmetric distributions. An example for L = 10 is shown in Fig. 11(e). 



4.4.2 Odd L: 



The lowest lying state has N+ = N- = and is the highest weight state of a [0, 1] multiplet. 
The sets {uj\j = 1, ... L — 1} and {7j|j = 1, ... L — 1} are closed under complex conjugation. 
Hence, after a suitable relabelling of roots we have 7^ = Uj, for j = 1, ... L — 1. Similarly as for 
the B = = state for even L this state can be considered as a collection of "2 — 2" narrow 
strings and the Bethe ansatz equations can be cast in the form (B.16) with N = L — 1. 
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Figure 11: Strange string solutions of the Bethe ansatz equations (4.6) in the octet sector for a system 
of size L = 10 with different AN = — N-: circles (triangles) denote uj (7a)- (a) AN = 1, (b) 
AN = 3, (c) AAT = 5 , (d) and (e) AN = 7, (f) N+ = 9,N-= 0. 

• A sequence of low energy [0, 1] octet highest weight states is obtained by taking N+ = + n, 
7V_ = ^ - n, n = 1,2, .... 

In the following section we carry out a finite-size scaling analysis for the states we have just described. 
B.3 Low-lying excitations in the thermodynamic Hmit 

In this Appendix we derive the coupled integral equations (4.17) from the algebraic equations (4.16). 
We start by taking the logarithm of (4.16) 

0{xj) + ^(|) = ^ + i E ^(^^) + 1 E ^(^) ' i = 1. • • • > 

fe=l fe=l 

(B.18) 

0{yj) + ^(|) = ^ + 1 E ^(^) + z E ^(^) ' = i> • • - 

fe=i fe=i 
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where 9{x) =2 arctan(2x) and Ij {Jj) are integer or half integer numbers characterizing the state 
uniquely. Following Yang and Yang we introduce counting functions z± {v) by 

z+{x) = e{x) + 6{-) -J 2^ ^(^^) ' 

k=l k=l 

(B.19) 

, N- N+ 

z. (X) = 0iy) + ^(|) - X E ^(^) - I E ^(^) ' 

k=l fe=l 

By definition, the counting function evaluated at a root of equations (B.18) yields 27r times the cor- 
responding integer numbers divided by L, e.g. z+{xj) = 2TrIj/L. Taking the derivative of (B.19) and 
then L — > oo we obtain the integral equations (4.17) for the densities of roots p± and holes p±. by using 
that 

2n[p1{x)+p±{x)] = ^^. (B.20) 
To proceed further, we restrict ourselves to states with symmetric distributions of the integers, e.g. 

Ij = l{-N+-l)+j, j = l,...,iV+ (B.21) 

and similarly for the Jfc. For these distributions we have 

p±{x) = 0, for \x\ < A± , (B.22) 

where A± are fixed by the conditions z±{A±) = N±/L. The integral equations for the root densities 
may then be simplified to 

rA+ rA- 
p+{x)+ dx'ttiix - x')p+{x') + 2 dy'a2{x-y')p-{y') = ai{x) + a3{x) , 

J-A+ J-A- 
fA- pA+ 

p-{y) + / dy'ai{y - y')p-{y') + 2 / dy'a2{y - x')p+{x') = a^{y) + as{y) , (B.23) 

J-A- J-A+ 

where 

1 2n 
an{u) = - 



TT 4^2 + 

The energy per site of such states in the thermodynamic limit is obtained from (4.6 

N. ^ JV_ 



1 1 1 

lim -En N_ = -27r lim < - ai(a;fe) + a3(a;fe) + - ai(yfe) + a3(yfe) 

— ►OO Li L—fOO h '—^ Li '—^ 

y fc=i k=\ 

fA„ 

= -27r^ / p<,(a;)[ai(u)+a3(w)]. (B.24) 



cr=± • 



B.4 Low- lying excitations for large, finite lattice sizes 
B.4.1 [0,-5, 5,0] states in the "sector" {7^} = {uj) for even L 

The considerations in the previous section are straightforwardly applied to the low-energy states with 
[B, J^] = [0, 0] in the indecomposable [0, — 5, ^, 0] multiplets for even L. For these we have to consider 
solutions satisfying 

7j=Uj, j = l,...,L, (B.25) 

which consist of degenerate type-± strange strings. In the thermodynamic limit the corresponding root 
distributions functions are given by (4.24) and hence we may use the counting function 

z{x) = 2 arctan (e'^^) - | . (B.26) 
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L 




L{El l - Eo) 


8 

10 

128 

256 

512 

1024 

2048 

4096 


-31.1255373481719104 

-39.3093986078541491 

-511.95141708922984 

-1023.9761543904651 

-2047.9882636363875 

-4095.9942108167784 

-8191.9971392992984 

-16383.998584343625 


6.9957012146247166 
6.9060139214585092 
6.2186125785810873 
6.1044760409276932 
6.0090181695995852 
5.9281236189417541 
5.8587150368839502 
5.7985285085172748 


oo 




4.95(2) (extrapolated) 



Tabic 7: Finite size results for the energy of the low-lying zero-momentum indecomposable consisting 
of two-strings with {uj} = {'Jj}. 



to obtain initial values (4.15) for a numerical solution of the Bethe ansatz equations (4.6): using 

Xj =yj = - IntauTT ( ^ + \ ] (B.27) 

TT \ L 4 J 

with 

Ij=-\{L + 2)+j, j = l,...,^ (B.28) 

as starting values a small number of iterations is sufficient to obtain a numerical solution of the Bcthc 
ansatz equations for lattices with more than 10.000 sites or L « 5000. We present our results on the 
finite size energies of an indecomposable in Table 7. The scaling behaviour as a function of 1/ log(L) 
is shown in Figure 5. 

B.4.2 Octet states in the "sector" {7^} = {uj} for odd L 

The same approach can be used to obtain the finite size energies of the lowest octet states of the odd 
L systems: these states are characterized by {L — l)/2 degenerate strange strings of each type -|- and 
— . For the numerical solution of the Bethe ansatz equations we use again (B.27) but with 

= ~{L + 1) + j, i = 1, . . . , ^ . (B.29) 

Table 8 summarizes results for the finite-size energies of the lowest-lying state in the octet sector. The 
momentum of this state is zero. For the extrapolation to L = 00 we assume the corrections of LAE to 
vanish as a function of l/log(i). The scaling of the energies together with the result of the rational 
function extrapolation is shown in Figure 5. 

B.4.3 More octet states built from strange strings 

Octet states from non-dcgcncratc strange string solutions exist for both even and odd system lengths 
and are described by Bethe ansatz equations with 2{L — 1) rapidities uj and arranged in A^_|_ {N-) 
strange strings of type -|- (— ) with + N- = L—1. We have solved the Bethe ansatz equations (4.6) 
for such states numerically for {AN = — N-): 

• L even: AN = 1,3,5,7 

• L odd: AN = 2 (AA'' = corresponds to the state studied in Section B.4.2) 

Starting values for the numerical solution of the Bethe ansatz equations for arc obtained from counting 
functions similar to (B.26) with proper symmetric choice of the quantum numbers Ij, Jk- The numerical 
finite size spectrum of these states is given in Table 9. The energy differences to the ground state as a 
function of 1/logL are plotted in Fig. 5. 
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Li 


77' 

nj L-1 L-i 

9. ■ 9. 


ijyil/L-i L-i — -lL/O) 


5 


-19.4230243905225 


2.88487804738772 


/ 


-Z r .ODDyDiDyoz4i>U 


o.UoizDoioozyyzy 


161 


-643.976352749708 


3.80720729705104 


321 


-1283.98786068577 


3.89671986852727 


641 


-2563.99379470064 


3.97759688996348 


1281 


-5123.99684487686 


4.04171273971125 


2561 


-10243.9983994774 


4.09893837823756 


5121 


-20483.9991898007 


4.14903019789338 


oo 




4.916(4) (extrapolated) 



Table 8: Finite size results for the energy of the low-lying zero-momentum octet consisting of two- 
strings with {uj} = {jj}. Rational function extrapolation has been used to obtain the i — > oo value. 



L 


Ell 1 

■> ' 9 


L(El l 1 — Eq) 


120 

240 

480 

960 

1920 

3840 


-479.96564417058130 
-959.98261219919220 
-1919.9912162303681 
-3839.9955678047386 
-7679.9977658894386 
-15359.998874847095 


4.1226995302440628 
4.1730721938711213 
4.2162094233208336 
4.2549074509588536 
4.2894922778941691 
4.3205871549434960 


oo 




4.91(2) (extrapolated) 


L 


El + 1 L-3 
■> > 9 


L(E L+l L-3 — £^0) 


121 

241 

481 

961 

1921 

3841 


-483.95681608386576 
-9.6397875255777501 
-1923.9895168083938 
-3843.9948148688195 
-7683.9974299991591 
-15363.998724015133 


5.2252538522429859 
5.1206335762225308 
5.0424151625950344 
4.9829110644573120 
4.9369716154324124 
4.9010578736833850 


oo 




4.937(3) (extrapolated) 


L 


Elj^-^l^2 


L{El^^ l_2 - Eg) 


128 

256 

512 

1024 

2048 

4096 


-511.94458966164291 
-1023.9737491964416 
-2047.9874366076488 
-4095.9939402783793 
-8191.9970594424680 
-16383.998566214848 


7.0925233097077580 
6.7202057109388988 
6.4324568838346750 
6.2051549395546317 
6.0222618254764871 
5.8727839789301150 


oo 




4.94(3) (extrapolated) 



Table 9: Finite size energies for the low-lying octets consisting of "strange strings" for L even and odd, 
respectively. 



42 



Based on our results we conjecture that for any fixed AN 

i!^^^ (^i(i-i+AJV),i(L-i-AJV) - Eo) « 4.93 , (B.30) 

i.e. we expect a macroscopic degeneracy of states in the finite size spectrum to order 1/L oi the super 
spin chain! 

B.5 The Fermi velocity 

In order to determine scaling dimensions we need to know the value of an important parameter, the 
Fermi velocity vp. Indeed, the Hamiltonian of our system is determined up to a scale only, and from 
conformal field theory wo only know, for instance, that the ground state energy Eq of the Hamiltonian 
on a lattice of L sites scales for L ^ oo like 

= eoo - + o(L ), (B.31) 

where Coo is the ground state energy density and c is the central charge of the CFT. The determination 
of Vp is straightforward is many systems solvable by Bethe ansatz, but not so here. We suggest that 
(we used this value in the main text to deal directly with a "properly normalized" Hamiltonian) 

w = TT. (B.32) 

Let us briefly list the evidence supporting (B.32). 

1. We can determine the dispersion of excited states over the lowest octet discussed above, that are 
described by the integer method (not all states are, as pointed out above). It follows from (4.8) 
and (4.9) that the bare energy and bare momentum of a 2-2 narrow string are given by 

eo{u) = -47r[ai(w)+a3(w)], (B.33) 
Po{u) = 2[e{u) + 9{u/3)], (B.34) 

where u denotes the real centre of the string. Given the integral equation (4.23) for the root 
density of the octet state, we may define a dressed energy e{u) and a dressed momentum p{u) as 
usual by 

/oo 
du' [2 a2{u-u') + aA{u-u')]e{u'), (B.35) 
-00 
1*00 

p{u) = Pq{u) — / du' [2 a2{u — u') + a4,{u — u')]p{u'). (B.36) 



These equations are easily solved by Fourier transform 

2tt 

e{u) = — — r- , p(u) = 2 arctan(sinh(7rM)). (B.37) 

cosh(7ru) 

We may eliminate the spectral parameter and obtain an expression of the dressed energy in terms 
of the dressed momentum directly 

e(p) = -27rcos(p/2). (B.38) 

The contribution of a 2-2 string with centre u to the ground state energy and momentum is then 
given by e(w) and p{u) respectively. Making a "hole" excitation by removing a 2-2 string with 
centre u from the ground state would change the energy and momentum by 

Eh = -e{u) , Ph = -p{u) , (B.39) 

and hence 

Eh{Ph) = 2TTCOs{Ph/2). (B.40) 
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Linearizing this dispersion around the "Fermi points" Ph = ikp = ±7r gives the Fermi velocity 

dEhiPh) 



V = 



dPh 



= TT. 



(B.41) 



2. Following [77] wc can determine v from functional equations for the largest eigenvalue A{u) of 
the transfer matrix by comparing them to crossing equations for relativistic S-matrices. From 
(4.11) with A = we obtain 

A{u) = A(l - u) , (B.42) 

which again implies that v = w. 

3. All of our finite-size spectra for small L ~ 5,7, 9, 11 systems are consistent with v = n as are the 

extrapolated values for very large systems. 

Now the scaling dimensions and conformal spins of the excitations considered above can be determined 
from the finite size spectra of the Hamitonian and the momentum operator 



x = h + h= = -^{Ec+4:L) 
L 

s = h-h = —P. 

2tt 



(B.43) 



B.6 "Antiperiodic" Boundary Conditions 

As discussed in the main text, "antiperiodic" boundary conditions play a crucial role in the identifica- 
tion of the continuum limit (even though they might not play any role in the physical network model, 
for instance). We have studied in some detail the sector M = N. The Hamiltonian is defined by taking 
the trace in (4.4) instead of the supertrace; observe that the boundary conditions break the s/(2/l) 
symmetry, but respect the bosonic spin-sM(2). The Bethe ansatz equations become 



Uj + 1 

Ui — i 



la +i 
la - i. 



-ll!^i^Ji±i, j = i,...,N 



N 



_YT la-u, + i ^ a = l,...,7V 



(B.44) 



The sector {uj} = {la} now reduces to the Takhtajan-Babujian chain with periodic boundary condi- 
tions. By diagonalizing the transfer matrix for small systems, we find that the ground state lies in the 
TB sector. For example for X = 3, 4, 5, 6 we have {ij = uj) 



L=3 


L=4 


L=5 


L=6 


E = -13.5335 


E = -17.4031 


E = -20.8890 


E = -24.8774 


ui = 0.3510 - 0.6258i 

U2 = —ul 
Us = 0.5219* 


Ul = 0.2958 - 0.5513* 

U2 = U* 

W3 = -0.2958- 0.5513* 

U4 = ^3 


Ul = 0.7564 + 0.8124i 
U3 = 0.1694-1-0.5980* 

U4 = — U3 

= -0.5013* 


Ul = 0.4385 + 0.5483* 
Us = -0.4385 + 0.5483* 

U4 = Us 

U5 = -0.5207i 

ue = ul 



The lowest-lying excited states above the TB ground state are formed by solutions of the Bethe ansatz 
equations (B.44) involving strange strings. In Fig. 6 we show the finite-size scaling of the TB ground 
state as well as two of these excited states. The numerical results indicate that both states become 
degenerate with the TB ground state in the L 00 limit. The scaling dimensions extrapolate to — ^ 
relative to the ground state with periodic boundary conditions. We conjecture that there again is an 
infinite degeneracy in the thermodynamic limit. 
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C Relation to the spin-1 Takhtajan-Babujian model 



Alcaraz and Martins [57] carried out a finite-size scaling analysis of the anisotropic spin-1 Takhtajan- 
Babujian (TB) chain based on the following equations 



[ sinh(7K+z]) |^ ,,,p(^^/Tf' sinh(7K-z^fe + i]) , . 

\ sinh(7K- -i])f '"P^*^^ n sinh(7K- - u, - i]) ' ^ ^^■^> 



(sin 27) 



2 L-n 

E 



2 cos 27 — cosh 27Uj 

L—n 

V = ^ 2arctan[tanh(7Uj)cot(7)] . (C.2) 
i=i 

The twist angle was taken to vary in the interval < $ < tt. Equations (C.1),(C.2) reduce to the 
equations of interest here if we take the limits $ — > tt, 7 — > and rescale £ by a factor of 4 and V by 
a factor of 2. The finite size spectrum in the TB chain in the limit $ — > tt, 7 ^ is of the form [57] 

£m,m' = Eo{^ = 0,L) + — + M + M') + o(L-i) , M, M' = 0, 1, 2, ... 

T^M,M' = Y{sa + M-M')+ o{L-') . (C.3) 

Here £'o($ = 0, L) is the energy of the ground state for zero twist-angle $ = and even L. It is equal 

to ^ 

£o($ = 0,L) = -L-^+o(L-i). (C.4) 
After the rescaling we then obtain the finite-size energies and momenta in the sZ(2/l) chain 



Em,m' = -4^L+ — {2xa-- +2M + 2M']+o{L-^) , M,M' ^0,1,2,... 



2^( _ 1 
L " 4 
27r 

Pm,m' = -J (2s« + 2M - 2M') + o{L-^) . (C.5) 

Here we have used that the Fermi velocity in the sl{2/\) chain is w = tt (B.32). This shows that the 
scaling dimensions in the TB model and the sZ(2/l) chain are related by 

• A doubling of the scaling dimensions in the s/(2/l) chain as compared to the TB model, combined 
with a shift by — j. 

• A doubling of the conformal spins, sectors). 

D Relation to the quantum transfer matrix of the supersym- 
metric t—J model 

The vertex model considered in this paper also appears in the quantum transfer matrix (QTM) approach 
to the thermodynamics of the integrable supersymmetric t-J chain: shifting the rapidities used in 
Ref. [54] (JKS) as 

1 I i , 

Vk o7/s + iu, Wj -Uj +IU+-, (D.l) 
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we obtain from (JKS.18) the BAE in our notation (4.6) with twisted boundary conditions 



when we set the free parameter A in our model to A = 4u+ 1. The twist entering (D.2) is related to the 
chemical potential fj, of the t-J model as exp(i^) = — exp(— /3/x) (/3 is the inverse temperature). Note 
that the QTM deals with the thermodynamic limit of the quantum chain, the system size N = 2L 
appearing in (D.2) is the so-called "Trotter number" appearing in the expression of the partition 
function by means of the QTM. The temperature dependence in the QTM approach enters via u = 



The eigenvalues of the QTM (.JKS.15/6) can bo identified with the eigenvalue (4.7) of the transfer 
matrix carrying the representation 3 in the matrix space 



provided that we choose the twist properly, i.e. exp(/3/i) = — 1. 

References 

[1] D. E. KhmeFnitskii, JETP Lett 38 (1983) 454; H. Levine, S. B. Libby and A. M. M. Pruisken, 
Nucl. Phys. B240 (1984) 30,49,71; A. M. M. Pruisken, Phys. Rev. Lett. 61 (1988) 1297; A. M. M. 

Pruisken, Nucl. Phys. B235 (1984) 277. 

[2] A. Altland and M. Zirnbauer, Phys. Rev. B55 (1997) 1142. 

[3] T. Senthil, J.B. Marston and M.P.A. Fisher, Phys. Rev. B60, 4245 (1999), cond-mat/9902062. 
[4] M. Janssen, M. Metzler and M. Zirnbauer, Phys. Rev. B59 (1999) 15836. 
[5] F. Evers, A. Mildenberger and A. D. Mirlin, Phys. Rev. B64 (2001) 241303. 
[6] L Affleck and D. Haldane, Phys. Rev. B36 (1987) 5291. 
[7] A. Zamolodchikov and Al. Zamolodchikov, Nucl. Phys. B379 (1992) 602. 
[8] K. B. Efetov, Adv. Phys. 32 (1983) 53. 
[9] H. A. Weidenmiiller, Nucl. Phys. B290 (1987) 87. 
[10] M. Zirnbauer, J. Math. Phys. 37 (1996) 4986. 

[11] M. Zirnbauer, "Conformal field theory of the integer quantum Hall plateau transition", hep- 
th/9905054. 

[12] N. Berkovits, C. Vafa and E. Witten, JHEP 9903 (1999) 018. 

[13] M. Bcrshadsky, S. Zhukov and A. Vaintrob, Nucl. Phys. B559 (1999) 205. 

[14] M. Zirnbauer, Ann. der Physik 3 (1994) 513; D. H. Lee, Phys. Rev. B50 (1994) 10788; J. Kondev 
and J. B. Marston, Nucl. Phys. B497 (1997) 639. 

[15] LA. Gruzberg, A.W.W. Ludwig and N. Read, Phys. Rev. Lett. 82, 4524 (1999), cond-mat/9902063. 

[16] N. Read and H. Saleur, Nucl. Phys. B613, 409 (2001). 




(D.2) 



-4/3/7V. 



^jKs{v) = -{v'^+ M^)^ A{2i{v + iu)) 



(D.3) 



46 



[17] R. Gade, J. Phys. A31 (1998) 4909. 

[18] I. Affleck, in "Fields, Strings and Critical Phenomena" , Les Houches 1988 
[19] I. Affleck, Nucl. Phys. B305 (1988) 582. 

[20] P. Goddard, D. Olive and G. Waterson, Comm. Math. Phys. 112 (1987) 591. 

[21] A. W.W. Ludwig, "A free field representation of the OSP{2/2) current algebra at level k = —2 
and Dirac fermions in a random SU{2) gauge potential", cond-mat/0012189 

[22] D. Bernard, "Conformal field theory applied to 2D disordered systems", hep-th/9509137; M.J. 
Bhaseen, J.S. Caux, I. Kogan and A. M. Tsvelik, Nucl. Phys. B618 (2001) 465. 

[23] N. Read and H. Saleur, to appear. 

[24] M. Marcu, .J. Math. Phys. 21 (1980) 1277;1284. 

[25] D. Leites, Acad. Nauk. CCCP (1983) 764 (in english). 

[26] J. Germoni, J. of Algebra 209 (1998) 367. 

[27] J. Germoni, "Representations indccomposables des superalgebres de Lie spcciales lineaires", Pub- 
lications de rinstitut de Recherche Mathcmatique Avancec, Universite de Strasbourg. 

[28] D. J. Benson, "Representations and cohomology", Cambridge Studies in Adv. Math. 30 

[29] H. G. Kausch, Nucl. Phys. B583 (2000) 513; M. Gaberdiel and H. G. Kausch, Nucl. Phys. B538 
(1999) 631. 

[30] C. de C. Chamon, C. Mudry and X.G. Wen, Phys. Rev. B53 (1996) R7638. 
[31] J. Jacobsen and H. Saleur, to appear. 

[32] F. Lesagc, P. Mathieu, J. Rasmussen and H. Saleur, Nucl. Phys. B647 (2002) 363. 

[33] M. Bocquet, D. Serban and M. Zirnbauer, Nucl. Phys. B578 (2000) 628. 

[34] P. Bowcock, M. Hayes and A. Taormina, Nucl. Phys. B510 (1998) 739. 

[35] P. di Francesco, P. Mathieu, Senechal, "Conformal field theory", Springer. 

[36] T. Eguchi and A. Taormina, Phys. Lett. B200 (1988) 315; Phys. Lett. B210 (1988) 125. 

[37] A.M. Semikhatov, A. Taormina and L Yu Tipunin, math.QA/0311314. 

[38] J. Links and A. Foerster, J. Phys. A 32, 147 (1999), cond-mat/9806129. 

[39] J. Abad and M. Rios, J. Phys. A 32, 3535 (1999), cond-mat/9806106. 

[40] R.M. Gade, J. Phys. A 32, 7071 (1999). 

[41] S. Derkachov, D. Karakhanyan and R. Kirschncr, Nucl. Phys. B583, 691 (2000). 

[42] M. Schemicrt, W. Nahm and V. Rittenberg, J. Math. Phys. 18, 146 (1977). 

[43] P. P. Kulish and E. K. Sklyanin, J. Sov. Math. 19, 1596 (1982). 

[44] P. P. Kulish, J. Sov. Math. 35 2648 (1985). 

[45] C.K. Lai, J. Math. Phys. 15, 1675 (1974); 

[46] B. Sutherland, Phys. Rev. B 12, 3795 (1975). 

[47] P. Schlottmann, Phys. Rev. B 36, 5177 (1987). 



47 



[48] F. H. L. Essler and V. E. Korepin, Phys. Rev. B 46, 9147 (1992). 
[49] A. Foerstcr and M. Karowski, Nucl. Phys. B396, 611 (1993). 

[50] F.H.L. Essler, V.E. Korepin, and K. Schoutens, Int. Jour. Mod. Phys. B8, 3205 (1994). 

[51] F.H.L. Essler and V.E. Korepin, Int. Jour. Mod. Phys. B8, 3243 (1994). 

[52] M. P. Pfannmiillcr and H. Frahm, Nucl. Phys. 479 [FS], 575 (1996); 
M. P. Pfannmiiller and H. Frahm, J. Phys. A30, L543 (1997); 
H. Frahm and M. P. Pfannmiiller, Zap. Nauch. Semin. POMI 251, 94 (1998). 

[53] H. Frahm, Nucl. Phys. 559 [FS], 613-636 (1999). 

[54] G. Jiittncr, A. KKimpcr and J. Suzuki, Nucl. Phys. B 487, 650 (1997), cond-mat/9611058. 

[55] H. Salcur and B. Wchcfritz-Kaufmann, Phys. Lett. B481, 419 (2000). 

[56] L. Takhtajan, Phys. Lett. A87, 479 (1982), 
H. Babujian, Nucl. Phys. B215, 317 (1983). 

[57] F.C. Alcaraz and M.J. Martins, J. Phys. A23, 1439 (1990), ibid A22, 1829 (1989). 

[58] H. Salcur and B. Wchcfritz-Kaufmann, Nucl. Phys. B663 (2003) 443. 

[59] P. di Francesco, H. Salcur and J.B. Zuber, Nucl. Phys. B300 (1988) 393. 

[60] J. Suzuki, J. Phys.A21, L1175 (1988). 

[61] H. de Vega and E. Lopes, Nucl. Phys. B362 (1991) 261. 

[62] H. Frahm and N.-C. Yu, J. Phys. A 23 (1990) 2115. 

[63] H. Frahm, N.-C. Yu and M. Fowler, Nucl. Phys. B 336 (1990) 396. 

[64] F.H.L. Essler, V.E. Korepin and K. Schoutens, J. Phys. A25, 4115 (1992). 

[65] F. Woynarovich, "On the symmetry of excitations in SU{2) Bethe ansatz systems", cond- 
mat/9812415. 

[66] P. Martin, "Potts models" , World Scientific 

[67] C.W. Curtis and I. Reiner, "Representation theory of finite groups and associative algebras", 

Wiley Interscience 1962. 

[68] O. Mathieu, in Adv. Studies in Pure Math., 26 (2000) 145. 

[69] V. Pasquier and H. Saleur, Nucl. Phys. B330 (1990) 523. 

[70] V. Chari and A. Pressley, "A guide to quantum groups" , Cambridge. 

[71] A. V. Razumov and Yu G. Stroganov, J. Phys. A34 (2001) 3185. 

[72] L.D. Faddeev and L. Takhtajan, J. Sov. Math. 24, (1984), 
A.N. Kirillov, J. Sov. Math. 34, 2298 (1985), 

F.H.L. Essler, V.E. Korepin and K. Schoutens, Phys. Rev. Lett. 67, 3848 (1991). 

[73] F.H.L. Essler, V.E. Korepin and K. Schoutens, Int. Jour. Mod. Phys. B8, 3205 (1994), 
K. Schoutens, Nucl. Phys. B413, 675 (1994). 

[74] G. Bediirftig and H. Frahm, J. Phys. A 28, 4453 (1995). 

[75] G. Bediirftig, F.H.L. Essler and H. Frahm, Phys. Rev. Lett. 77, 5098 (1996), 
A. Foerster, J. Links and A.P. Tonel, Nucl. Phys. B552, 707 (1999). 



48 



[76] H. Frahm, M. P. Pfannmiillcr and A. M. Tsvclik, Phys. Rev. Lett. 81, 2116 (1998). 
[77] M.J. Martins, B. Nienhuis and R. Rietman, Phys. Rev. Lett. 81, 504 (1998). 
[78] H. Saleur, Nucl. Phys. B578, 552 (2000). 



49 



